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Introduction 



In the French academic system, the Habilitation a diriger des recherches, or HDR, is a title that 
allows one to supervise the research of PhD students. Therefore, to obtain it, the candidate is 
required to prove that he has a wide enough view of his research subjects, and the capability to build 
a comprehensive research project, such that he can guide his research, and that of others, in the short 
and long runs. At Ecole Normale Superieure, the length of the thesis is supposed to be around 40 
pages. 

Given these constraints, I decided to focus on one of my research subjects only, on which I have 
been working continuously over the last eight years, and where I strongly contributed to the design of 
the research. This is the theory of the glass and jamming transitions of hard spheres. On this topic, 
I recently wrote a very detailed review paper [1] as well as a long research paper [2]. I refer to these 
two papers for most of the technical details. Here I only tried to explain this subject in a concise and 
self-contained way, highlighting what I consider the main motivations for its study and trying to put 
my work in a broader perspective. Basically I tried to tell this story in the same way as I would tell 
it to a student who contacted me to start a PhD thesis on these topics. 

The glass problem has a reputation of being controversial. It has even been said that "There 
are more theories of the glass transition than there are theorists who propose them" [3]. Part of 
this controversy is fully justified. The problem is indeed extremely challenging, and it has long 
been difficult for theory to provide precise predictions that could be experimentally tested, in order to 
obtain conclusive answers [4] . A lot of progress has nonetheless been achieved thanks to the impressive 
numerical and experimental developments of the last couple of decades, and theory is rapidly evolving 
in the aim to match this effort. 

Yet part of the controversy may also be philosophical 1 . There seems to be some disagreement 
about what one expects from a theory of the glass transition. Before proceeding, I thus find it useful 
to briefly review this issue, in order to clarify my scientific intent. I will use a convenient definition 
of "scientific theory" that has been discussed by Lucio Russo in [5, Sec 1.3] (together with a nice 
discussion of its limits that is not reproduced here). According to this definition, a "scientific theory" 
must have the three following properties: 

1. Its statements do not concern concrete objects pertaining to the real world, but specific abstract 
objects. 

2. It has a deductive structure: it is made by a few postulates concerning its objects, and by a 
method to derive from them a potentially infinite number of consequences. 

3. Its application to the real world is based on a series of "correspondence rules" between the 
objects of the theory and those of the real world. 

According to Russo, a useful criterion to determine whether a theory has these properties is to check 
if one can compile a collection of exercises, that can be solved within the theory. Indeed, solving a 

X I also fear that another part of the problem may be systemic. The institutional organization of scientific research 
seems sometimes designed to promote the radicalization of scientific discussions. Most of us experience a strong pressure 
to highlight the difference between one's theory and those of others, in order to publish our papers, raise more funds, etc. 
For example, a journalist was recently interested in writing about the publication of one of my papers for a prestigious 
journal, but the journal's editor refused to publish the news piece when he was told that our results did not generate 
sufficient controversy. Our results partially agreed with an existing theory, so they were considered not interesting 
enough! 
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problem in the context of the theory is nothing else but an (arbitrarily difficult) "exercise". This 
definition is of course very restrictive, as it rules out all phenomenological and empirical knowledge 
in science. The aim here is not to claim any superiority of "scientific theories" with respect to other 
more empirical or descriptive forms of science, but to provide a precise definition of a theory, in order 
to avoid useless controversies based on disagreement on what kind of "theory" one is looking for. 
"Scientific theories" are powerful for two reasons: 

(i) Working on two parallel levels (the "model" and the real world) allows for a very flexible way of 
reasoning, and in particular one can guarantee the "truth" of scientific statements by limiting 
them to the domain of the theory. 

(ii) The theory can be extended, by using the deductive method and introducing new correspondence 
rules, to treat situations that were not a priori included in the initial objectives for which the 
theory was developed. 

At the same time, it is important to keep in mind that any "scientific theory" has a limited utility, as 
in general it can only be used to model phenomena that are not too "far" from those that motivated 
its elaboration. Theories that become inadequate to describe a new phenomenology must be, for this 
reason, substituted; they remain however, according to our definition, "scientific theories" and one can 
continue to use them in their domain of validity 2 [5] . This last statement is particularly important 
because it reminds us that the theory of a given class of phenomena -e.g. the glass transition- will 
never exist 3 : scientific theories are never unique nor everlasting. They are models of reality, and there 
is no problem in using different models of the same phenomenon, and in replacing current models 
with more powerful ones when they are found. 

I will argue in this thesis that the so-called Random First Order Transition, or RFOT, theory -on 
which my work is based- is indeed a scientific theory of the glass transition according to the definition 
above. In particular, I will present a striking example of the usefulness of property (ii) above: RFOT 
theory, originally born as a theory of the thermal glass transition, can be surprisingly extended to 
describe the physics of the athermal jamming transition that happens inside the glass phase. One of 
the main criticism to RFOT theory, which has been a source of endless debates, is that it predicts 
the existence of a phase transition (the Kauzmann transition) that, in the real world, is at best 
unobservable and at worst non existent. However, according to the definition above, this is not a big 
issue: it just means that this "abstract object" is not related by a correspondence rule to a real world 
phenomenon, and is just a tool that is used within the theory to model other, observable, phenomena 4 . 
Other theories of the glass transition, for instance those based on the idea of dynamic facilitation, 
provide in my opinion nice complementary views on the problem. Overall, exciting progresses have 
been made, and even if we don't have -and will never have- the theory of the glass transition, we 
learned how to build consistent models of many phenomena related to it and there is reason to hope 
that even bigger progresses will be made in the near future. 

In the following, Chap. 1 provides an introductory account of RFOT theory, while Chap. 2 reviews 
my original work on this problem. A final comment on the title of this thesis is needed here. The 
approach to the glass transition that I will present is really similar in spirit to the one of the celebrated 
book Theory of simple liquids by Hansen and McDonald [6] . This approach is in fact based on a natural 
extension of the techniques used in that book to describe "simple" liquids -monoatomic ones, modeled 
as hard spheres or point particles interacting through simple potentials such as the Lennard-Jones- to 
the glassy states that are obtained by cooling or compressing such liquids. I hope that soon enough 
the theory will be so much developed that someone will write a book with this title, which would be 
the natural "volume 2" of a series on the amorphous states of matter, "volume 1" being [6]. 



2 For instance, Newtonian mechanics did not become useless once quantum mechanics was developed, and the fact 
that it gives incorrect predictions -e.g. the instability of the hydrogen atom- does not mean that it is plain wrong. 

^Contemporary popular literature on science sometimes forgets this point, e.g. when talking about particle physics. 

4 A similar situation was present in the standard model of particle physics until a few months ago, when the discovery 
of the Higgs boson was announced. But even if the Higgs boson was not discovered, the standard model would have 
retained its validity as a theory of low-energy particle physics. The Higgs boson would have been just an object of the 
theory, without any real world correspondent. 



Chapter 1 

A theory of the glass transition 



1.1 Basic facts 

1.1.1 Viscosity and structural relaxation time 

Although liquids normally crystallize on cooling, there are members of all liquids types (including 
molecular, polymeric, ionic and metallic) that can be supercooled below the melting temperature T m 
and then solidify at some temperature T gi the glass transition temperature. The viscosity rj(T) of 
the liquid increases continuously but very fast below T m and at some point reaches values so high 
that the liquid does not flow anymore and can be considered a solid for all practical purposes: at 
low temperatures, an amorphous solid phase is observed. The temperature T g marking the transition 
between the liquid and the glass is often conventionally defined by the condition f](T g ) = 10 13 Poise, 
but many other definitions are possible. As an example of this phenomenon, in Fig. 1.1 the viscosity 
of many glass forming liquids is reported as a function of the temperature. Following Angell [7], the 

quantity log 10 is reported as a function of T g /T. The viscosity increases of about 17 orders of 

magnitude on decreasing the temperature by a factor 2. Note that because the increase of viscosity is 
so fast, the dependence of T g on the particular value of viscosity (10 13 Poise) which is chosen to define 
it is very weak. The viscosity around T g is often described by the Vogel-Fulcher-Tamman (VFT) law, 

V(X) = fee 7 ^ , (1.1) 

where f^, A and To are system-dependent parameters. If To = this relation reduces to the Arrhenius 
law; otherwise, the extrapolation of the viscosity below T g leads to a divergence at T — Tq. 

The viscosity is related to the structural relaxation time T a by the Maxwell relation, n = GooT" a , 
where Goo is the infinite-frequency shear modulus of the liquid. The structural relaxation time is 
related to the decorrelation of density fluctuations. In glass forming liquids, for T m S> T > T g , the 
decorrelation of density fluctuations happens on two well separated time scales: a "fast" time scale 
(to <~ 10~ 12 s), which is related to vibrations of the particles around the disordered instantaneous 
positions, and a "slow" time scale r Q , which is related to cooperative rearrangements of the disordered 
structure around which the fast vibrations take place. Through the Maxwell relation, the fast increase 
of viscosity around T g is then related to a marked slowing down of the structural dynamics; usually, 
at T g one has r Q ~ 100s ~ 10 14 r , while in the liquid phase r Q ~ fis. The structural relaxation time, 
obtained from dielectric relaxation data, of some fragile glass forming liquids is reported in the right 
panel of Fig. 1.1. The behavior of r Q (T) is also described by a VFT law with an apparent divergence 
at T = To. This leads to the interpretation of T as a temperature at which a structural arrest takes 
place. 

Whether such a divergence really exist in glasses is hotly debated. The reason is that because 
of the strong divergence of Eq. (1.1), the T g below which the system cannot be equilibrated is quite 
far from T : typically, T g rts |T . Hence it is very hard to approach the putative divergence point 
while being at equilibrium. For comparison, in standard critical phenomena the relaxation time goes 
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Figure 1.1: (Left, from [7]) Viscosity data for many glass forming liquids. The logarithm of the 
viscosity measured in Poise is reported as a function of T g /T. The (calorimetric) T g is defined as the 
temperature at which the enthalpy relaxation time is ~ 200s, and its value is reported in parenthesis 
in the key of the figure. Note that for some systems the value of the calorimetric T g does not satisfy 
exactly the condition f](T g ) = 10 13 Poise. Fragility is the slope of the curves in T g /T = 1. (Right, 
from [8]) Structural relaxation time obtained from dielectric relaxation measurements. The dashed 
line indicates Arrhenius behavior. The values of T g , obtained from r a {T g ) — 100s, are reported in the 
key. 



as t~ t |T/T c — l\~ zl/ with an exponent zv <~ 1 [9]. This means that if one is able to equilibrate on 
a time scale r/r = 10 14 , then one can approach the critical point at a distance \T — T c \ <~ 10~ 14 T C . 
Hence, in this case one can equilibrate the system arbitrarily close to the critical point, the only 
limitation being the precision on temperature control. 



1.1.2 Correlation functions 

Consider a simple glass former made by N particles with positions xj. An important observation is 
that, if pk = N~ x / 2 e lkXi is the Fourier transform of the density profile, static correlation functions 
such as the structure factor 

S(k) = lim I = lim (p k p_ k ) (1.2) 

N->oo iv *■ — ' N— too 

jl 

do not exhibit any particular change of behavior at the glass transition: they are very smooth functions 
of temperature in all the range of temperatures above and below T g . 

The most spectacular signature of the incoming transition can be found in time-dependent corre- 
lations, such as the mean-square displacement 



(1.3) 
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Figure 1.2: Dynamical correlations from a numerical simulation of a binary mixture of particles inter- 
acting through Lennard- Jones potentials, with parameters chosen in such a way to avoid crystallization 
and obtain a good glass former. (Left, from [11]) Incoherent scattering function at a fixed wavevector 
(here denoted by q) as a function of time and several temperatures above the glass transition. (Right, 
from [10]) The mean square displacement as a function of time for the same temperatures. The plateau 
at intermediate times is seen to appear at the lower temperatures. 

and the so-called coherent and incoherent scattering functions, respectively defined by: 
F{k,t) = lim l£yM*,(t)-*,(0)]) = Um (p k ( t )p_ k (0)) , 

! * (1-4) 
F s (k,t)= lim — J2{e iklXj{t) ~ Xj{0)] ) ■ 



j 

Here the average is over the dynamical process, at equilibrium in the liquid phase. 

An example of such correlations is given in Fig. 1.2, for a numerical simulation of a Lennard- 
Jones system [10, 11]. Close to the glass transition, time-dependent correlations display two distinct 
relaxations, separated by a region where they depend weakly on time (a plateau). Based on these 
results, a common pictorial interpretation of the dynamics of glass forming liquids above T g is the 
following: for short times the particles are "caged" by their neighbors and vibrate around a local 
structure on a nanometric scale; the structural relaxation is then interpreted as a slow cooperative 
rearrangement of the cages. Note that on the time scale of the structural relaxation time T a , the root 
mean square displacement of the particles is of the order of 10% of the particle radius, so one cannot 
think to the structural relaxation as a process of single-particle "jumps" between adjacent cages. 

At T g and below, the structural relaxation time is much larger than the experimentally or numer- 
ically accessible window and the plateau extends to the largest accessible times. Hence one can define 
a caging order parameter as 

A= lim (Ar 2 (t)) , 
f(k)= lim F(k,t)/S(k) , 

t— >oo 

where the limit t — > oo has to be intended as taking the largest accessible times. A is often called cage 
radius while f(k) is often called non-ergodicity parameter. If the relaxation time were really to diverge 
at some To, then this ideal glass transition would be characterized by the (discontinuous) appearance 
of truly finite long-time limits of dynamical correlations, signaling the complete structural arrest. 

It is also interesting to consider the out of equilibrium dynamics following a fast quench into the 
glass phase. Suppose that an equilibrium liquid prepared at some temperature Tj > T g is rapidly cooled 
(quenched) at a temperature Tf < T g . Define t = as the time at which the cooling procedure ends 
and T — Tf. One can wait a time t w after the quench and then measure the dynamical correlations 
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F(k;t,t w ) — liniAr^oo (p^ (t + t w )p-k(t w )). Then, it is observed that increasing t w (the age of the 
system) has roughly the same effect as decreasing temperature at equilibrium. One observes a two- 
step relaxation with a structural relaxation time that increases with t w until the structure becomes 
arrested. This process is known as aging and it is another important characterization of the glass 
phase. 

1.1.3 Configurational entropy 

The idea that the dynamics in the supercooled phase is separated in a fast intra-cage motion and in 
a slow cooperative rearrangement of the structure suggests to split the total entropy of the liquid in a 
"vibrational" contribution, related to the volume of the cages, and a "configurational" contribution, 
that counts the number of different disordered structures that the liquid can assume [12]: 

S lig (T) ~ S vib (T) + S C (T) . (1.6) 

To estimate the vibrational contribution to the entropy of the liquid, it is often assumed that it is 
roughly of the order of the entropy of the corresponding crystal. Despite the fact that this idea is plain 
wrong for hard spheres and similar system where excluded volume effects play an important role, this 
approximation is not so bad in systems -such as the Lennard- Jones potential- where interactions are 
smoother and have longer range. It is then possible to estimate the configurational entropy S C (T) as 

/T m jrpl 
[Q iq (T>) - C cryst (T')} , (1.7) 

where AS m = Su q (T m ) — S crys t(T m ) is the entropy difference between the liquid and the crystal at the 
melting temperature T m , and C(T) = is the specific heat. Note that in experiments one usually 
works at constant pressure, C = C p , while in numerical simulations and in theoretical computations 
one usually works at constant volume, C — C v . 

In Fig. 1.3 the estimate of S c , obtained from calorimetric measurements of the specific heat and 
using Eq. (1.7), is reported for four different fragile glass formers. Below T g the liquid falls out of 
equilibrium as the structural relaxation time becomes of the order of the experimental time scale 
(~ 100s). This means that the structural rearrangements are "frozen" on the experimental time scale 
and the only contribution to the specific heat comes from the intra-cage vibrational motion; in this 
situation the specific heat of the liquid becomes of the order of the crystal's one and S C (T) approaches 
a constant value. However, one can ask what would happen if the time scale of the experiment were 
much bigger, say 10 6 s. In this case, the glass transition temperature T g would be lower and the 
plateau would be reached at smaller values of S c . If one assumes to be able to perform an infinitely 
slow experiment, one can imagine to follow the extrapolation of the data collected above T g to lower 
temperatures. For fragile liquids, it is found that a good extrapolation is 

Sc(T) = S 00 (l-?f) , (1.8) 

where the parameters Soo and Tk are fitted from the data above T g . This extrapolation is reported 
as a full line in Fig. 1.3. 

The outcome of this procedure is that the configurational entropy seems to vanish at a finite 
temperature Tk- Because S c counts the number of different structures that the liquid can access, it is 
not expected to become negative. A possible explanation of this paradoxical behavior was proposed 
by Kauzmann [12], who argued that at some temperature between T g and Tk the free energy barrier 
for crystal nucleation becomes of the order of the free energy barrier between different structures of 
the liquid. This means that the time scale for crystal nucleation becomes of the order of the structural 
relaxation time r a of the liquid, and one cannot think anymore of an "equilibrium" liquid because 
crystallization will occur on the same time scale needed to equilibrate the liquid. The extrapolation 
of S C (T) down to Tk is then meaningless, and the paradox is solved. This argument has been recently 
reconsidered, see e.g. [14], and its implications are still under investigation. 

Alternatively, one can assume that the existence of the crystal is irrelevant, because crystallization 
can be in some way strongly inhibited: for instance, by considering binary mixtures [10, 11], or -in 
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T/K T/K T/K T/K 

Figure 1.3: (From [8] and references therein) Configurational entropy S C (T) of four fragile glass 
formers. The black squares are obtained from calorimctric measurements of the specific heat of the 
liquid and of the crystal, see Eq. (1.7). Below T g (reported in the key) the liquid falls out of equilibrium. 
The black line is the extrapolation according to Eq. (1.8) of the equilibrium data for T >T g below T g , 
that goes to zero at T — Tk- The open white circles are another estimate derived from the dielectric 
relaxation data of Fig. 1.1 using the Adam-Gibbs relation [13]. 

numerical simulations- by adding a potential term to the Hamiltonian that inhibits nucleation [15]. 
If crystallization is neglected, the extrapolation of S c suggests that at Tk a phase transition happens, 
at which the number of structures available to the liquid is no more exponential, as S c = 0, and 
the system is frozen in one amorphous structure which can be called an ideal glass. Below T K , the 
only contribution to the entropy of the ideal glass is the vibrational one, so the specific heat has a 
downward jump at Tk, corresponding to the freezing of some degrees of freedom at the transition. 
The transition is expected to be of second order from a thermodynamical point of view. 

An evidence supporting this picture is the fact that in almost all the fragile glass formers it is 
found that Tk ~ To. For instance, in [16] some 30 cases where T = Tk with an error of order 3% 
are reported (but the reliability of these extrapolation has been questioned in [17]). The equality of 
(non-zero) To & n d Tk implies that both the structural relaxation time and the viscosity diverge at 
Tk, so that the structures that are reached at T K are thermodynamically stable, being associated 
to an infinite structural relaxation time. Of course, this ideal glass transition, that would occur in 
equilibrium, would not be observable: at some temperature T g > Tk where T a (T g ) = T exp a real glass 
transition, freezing the system in a nonequilibrium amorphous state (a real glass), happens. The value 
of T g , as well as the properties of the nonequilibrium glass (density, structure, etc.) depend on the 
value of r exp , which is usually ~ 100s as already discussed. 

1.1.4 Dynamical heterogeneities 

An important question is whether the slow relaxation observed on approaching T g is produced by 
independent relaxation events that are due to locally high potential energy barrier, or by cooperative 
effects that make relaxation difficult. To quantify this, it is convenient to introduce a correlation 
function that measures how much structural relaxation is correlated in space. 

Consider the real-space density profiles p(x,t) = J2iLi $i x ~ x i{t)) at ti me zer0 an d at time t, re- 
spectively given by p(x, 0) and p(x, t). One can define a local similarity measure of these configurations 
as 

C(r,t)= Jdxf(x)p(r+^,t)p(r-^,0)-p 2 , (1.9) 

where f(x) is an arbitrary "smoothing" function of the density field with some short range tq. In 
experiments, f(x) could describe the resolution of the detection system and can be for instance a 
Gaussian of width r . Call C(t) = V^ 1 J dr(C(r,t)) the spatially and thermally averaged correlation 
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Figure 1.4: (From the cover of PNAS, September 8, 2009) Dynamical heterogeneity in a structural 
glass-forming liquid. Shown is a spacetime rendering of the equilibrium dynamics of a two-dimensional 
super-cooled fluid mixture of 10000 particles after a fraction of a structural relaxation time. Particles 
are colored according to their overlap with their initial positions: a particle that is displaced by more 
than one particle diameter is dark red; a particle that has no displacement is dark blue; intermediate 
colors in the visible spectrum coincide with intermediate displacements. The color variation illustrates 
significant dynamic heterogeneity. Yet the spatial arrangement of particles at a given time seems 
perfectly homogeneous when coarse grained over only one of two particle diameters. The juxtaposition 
shows that the dynamics of this system are highly correlated, but the structure is seemingly not. 

function. On approaching the glass transition, this function behaves exactly as the coherent scattering 
function defined in Eq. (1.4). 

The correlation function C(r, t) probes the dynamics in the vicinity of point r. Spatial correlations 
in the dynamics can be quantified by introducing the correlation function of C(r,t), i.e. a four-point 
dynamical correlation 

G 4 (r, t) = (C(r, t)C(0, t)) - (C(r, t)) ((7(0, t)) . (1.10) 

The latter decays as Gi(r : t) ~ exp(— r/£(i)) thus defining a "dynamical correlation length" £(t). In 
a variety of glass formers, it is found that grows when t approaches values corresponding to the 
plateau, and reaches its maximum £ on the scale of the structural relaxation. The value of £ is found 
to increase upon decreasing temperature towards T g [18]. 

A large value of £ indicates that if a given region in space is "mobile" -i.e. it has a structural 
relaxation faster than the average- then neighboring regions will also be mobile, and similarly for 
slow regions - those with a structural relaxation slower than the average. This fact indicates that 
dynamics is cooperative and characterized by a "facilitation" mechanism [19, 20, 21, 22]: a mobile 
region can speed up the dynamics of neighboring regions thanks to dynamical correlations. Often 
such cooperatively rearranging regions, whose presence is now well established, are called "dynamical 
heterogeneities" [18]. More pictorial measures of dynamical heterogeneities are possible, one of them 
is reported in Fig. 1.4. 

An important question is whether dynamical heterogeneities are purely dynamical in nature or 
if they are induced, or accompanied, by the growth of some static correlation length. As previously 
discussed, simple measures of structural correlations do not show any particular anomaly at the glass 
transition, while the dynamics displays a strong slowing down and increasing heterogeneity. There 
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t=40000 MCsteps t=1 000000 MCsteps 




Figure 1.5: (From [28]) Fluctuations of the overlap field, SC(r, t) = C(r, t) — C(t), for a two dimensional 
slice of the three-dimensional simulated system. Upper panels: unconstrained system. Lower panels: 
constrained system (C = 0.25). Left panels: t = r a . Right panels: large times. Phase separation is 
observed in the constrained system. 

have been several attempts to define and measure such a static correlation length, leading in particular 
to the introduction of the so-called point-to-set correlation length [23, 24, 25] which has been measured 
e.g. in [26, 27]. Perhaps the most spectacular result is the one obtained in [28]. Recall that C(t) 
measures the average "similarity" (overlap) between the density profile at time t = and t. In this 
numerical simulation of a Lennard- Jones-like glass former, a constraint has been added to the system, 
imposing that C{t) > C, a fixed constant. In other word, the overlap with the initial configuration 
(which is just an equilibrium reference configuration) must be bigger than some prefixed constant. 
Hence, this constraint is static in nature: the system can visit all configurations with overlap larger 
than C according to its equilibrium distribution. The result of this investigation, reproduced in 
Fig. 1.5, show that the system will phase-separate into regions where the overlap is very low and 
regions where overlap is very high. This indicates that some static "surface tension" between regions 
of high and low mobility exists, and is responsible for dynamical heterogeneities. Similar results have 
been obtained by imposing the overlap constraint in different ways, see e.g. [29, 30]. 

1.1.5 Glass and jamming transitions 

Besides the glass transition, another "rigidity" transition has attracted a lot of attention: the so-called 
"jamming" transition of granular matter. Observing a jamming transition is an everyday experience, 
similarly to the glass transition. An athermal amorphous assembly of hard objects -such as nuts, 
oranges, tennis balls- will be mechanically stable, meaning that it can support finite stresses, if its 
density is large enough. In absence of friction, for three-dimensional spheres, typically around ~ 64% 
of space is filled, while in presence of friction the density can be lower. For comparison, in three 
dimensions the closest packing of an assembly of hard spheres is such that ~ 74% of space is occupied, 
corresponding to a bcc/fcc crystalline structure. 

Because granular matter is typically made by hard objects, we need to introduce a model system 
that can interpolate between finite temperature glasses and hard spheres. A very convenient one is the 
model of frictionless harmonic spheres: N spherical particles of diameter a are enclosed in a volume 
V in d spatial dimensions, and interact with a soft harmonic repulsion of finite range: 

v(r) = e(l-r/a) 2 6{o- ~r) , (1.11) 

where r > is the interparticle distance, 0(r) is the Heaviside step function and e controls the 
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Figure 1.6: A schematic phase diagram for the glassy states of the model defined by Eq. (1.11) showing 
both the glass transition line at T g (ip) and the jamming point at T = and <p = cpj. Below ipj, particles 
do not overlap, as in hard sphere configurations. Above ipj, there is a finite density of overlaps, and 
the energy and pressure are finite. In this region, hard sphere configurations cannot be found. 

strength of the repulsion. In this model, there is no force between particles if they do not overlap, 
while there is a harmonic repulsion if the particles overlap. This is supposed to mimic rigid particles 
with a finite elasticity. The model (1.11), originally proposed to describe wet foams, has become a 
paradigm in numerical studies of the T = jamming transition [31, 32, 33]. It has also been studied 
at finite temperatures [34, 35, 36], and finds experimental realizations in emulsions, soft colloids and 
grains [37]. The model has the two needed control parameters to explore the jamming phase diagram: 
the temperature T (expressed in units of e), and the fraction of volume occupied by the particles in the 
absence of overlap, the packing fraction if = pVd(o~/2), where V^(r) is the volume of a d-dimensional 
sphere of radius r and p = N/V is the density. 

Over the last decade, a large number of numerical observations have been reported for this 
model [32, 33]. A jamming transition is observed at T = at some critical volume fraction (pj, 
the density above which the packings carry a finite density of particle overlaps. This transition is pic- 
torially represented in Fig. 1.6. Numerically, the zero-temperature energy density, ecs, and pressure, 
P, are found to increase continuously from zero above ipj as power laws [31]. The pair correlation func- 
tion of the density fluctuations [6], g(r), develops singularities near <pj [38, 39], which are smoothed 
by thermal excitations [36]. In particular, g(l) = oo at ip = <pj and T = 0. This behavior implies 
that the density of contacts between particles, z, jumps discontinuously from to a finite value, Zj, at 
ifj, and increases further algebraically above (fij [31, 32]. Thus, the jamming transition appears as a 
phase transition taking place in the absence of thermal motion, with a very peculiar critical behavior 
and physical consequences observable experimentally in a large number of materials. 

To make a connection with the glass phenomenology discussed above, note that the liquid un- 
dergoes a glass transition at some temperature T g (<p). The temperature T g (ip) vanishes at a volume 
fraction ip g ~ 58% which corresponds to the glass transition of hard spheres. This is because when 
e>T, the harmonic spheres model reduces to hard spheres, provided the density is low enough that 
overlaps can be avoided. The point <p g is not the jamming transition: indeed, the ground state energy 
and pressure remain zero across ip g . Above ip g , the system enters at zero temperature a hard sphere 
glassy state. In this state, particles vibrate near well-defined (but random) positions, and the system 
is not yet jammed. Jamming happens when the glass reaches its close packing density, which is then 
identified with the jamming transition at (pj. In Fig. 1.6 the density interval [ip g ,ipj] ~ [0.58,0.64] is 
simply the amorphous analog of the interval tp € [(p m ,<p cp \ ~ [0.54,0.74] for ordered states of hard 
spheres where a compressible crystalline structure exists at thermal equilibrium. 

The crucial difference between the hard sphere glass phase for ip <E [cp g , ipj] and the jammed phase 
for ip > ifj is the following. At any finite temperature, both phases are rigid, in the sense that they 
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have finite elastic moduli, pressure, etc. However, the rigidity of the glass phase is entropic in nature: 
it is due to thermal vibrations of the particles, that induce collisions between neighboring particles 
and allow them to exchange momentum. Indeed, pressure and all elastic moduli are proportional to 
temperature in this regime, and they vanish at zero temperature [6, 40]. On the contrary, the rigidity 
of the jammed phase is due to overlaps between particles, that produce interaction forces even at 
strictly zero temperature. Above ipj, pressure and clastic moduli are finite even at T = [31, 32, 33]. 
In summary, ipj marks a transition inside the glass phase from entropic to mechanical rigidity. Below 
<Pj, some noise (e.g. due to thermal agitation) is needed to stabilize the system, that would otherwise 
be mechanically unstable. Above ipj, the system is rigid in absence of any noise. 

The main theoretical challenge in studying the jamming transition at ipj is that it happens deep 
inside the glass phase, hence out of equilibrium. Therefore, an accurate theoretical description of 
the glass phase is needed. However, once this is achieved, the jamming transition provides a unique 
opportunity to perform very detailed tests of theories of the glass structure. Indeed, while the structural 
signatures of the glass transition are quite elusive (no signature at all in the pair correlation g(r), weak 
singularities in the thermodynamics), the structural signatures of the jamming transition are extremely 
easy to detect numerically: at ipj marked changes happen in g(r), and pressure and energy have strong 
singularities. 



1.1.6 Summary 

In summary, the liquid-glass transition is a very strange transition, much different from standard 
critical phenomena. It is characterized by: 

• A strong increase of relaxation time upon lowering temperature gives the impression of a diver- 
gence at a finite temperature To. However the divergence is so strong that the system falls out 
of equilibrium at a T g which is quite far from T , in such a way that approaching the putative 
critical point while being at equilibrium is impossible. 

• A separation of time scales between fast vibrational dynamics (to ~ 10 _12 s) and slow structural 
relaxation r Q 3> tq. Dynamical correlations exhibit a characteristic plateau for tq <C t <C T a . 
At T g the length of the plateau exceeds the accessible time scales and the system is effectively 
frozen (in this sense, T g effectively acts as T , except for the fact that the system falls out of 
equilibrium). The value of mean square displacement at this plateau defines the cage radius A 
at the glass transition, while the long time limit of the coherent scattering function defines the 
non-ergodic parameter f(k), Eq. (1.5). 

• A strongly heterogeneous dynamics, with well-separated regions containing mobile or slow parti- 
cles. Although everybody agrees on the existence of these regions as defined from the dynamics, 
some have argued that they have a static origin, and have given correspondingly numerical 
evidence for a growing static correlation length. 

• The equilibrium configurational entropy of the liquid, S C (T), is found to decrease quickly around 
T g , and to be correlated with the increasing of relaxation time (as expressed by the Adam- 
Gibbs relation [13]). If one believes in extrapolations, the extrapolation of S C (T) to zero defines 
a temperature Tk which is found, for several systems, to be quite close to the extrapolated 
divergence of relaxation time at To. Moreover, the freezing out of structural degrees of freedom 
at T g induces a sudden drop in specific heat (or in compressibility for hard spheres). 

In addition, for systems with a finite radius of interaction, a jamming transition happens inside the 
glassy region of the temperature-density phase diagram, at which the rigidity of the glass changes in 
origin. A correct description of the glass phase should allow one to describe this transition and its 
geometrical and structural signatures. 
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Theory of second order PT (gas-liquid) 


RFOT theory of the liquid-glass transition 


Qualitative MFT 


Landau theory (1937) 
Spontaneous Z 2 symmetry breaking 
Scalar order parameter 
Critical slowing down 


Replica theory (Parisi, 1979; KTW, 1987) 
Spontaneous replica symmetry breaking 
Order parameter: overlap matrix q a b 
Dynamical transition "a la MCT" 


Quantitative MFT 
(exact for d — » oo) 


Liquid-gas: /3p/p = 1/(1 — pb) — j3ap 
(Van der Waals 1873) 
Magnetic: m = tanh(/3 Jm) 

(Curie- Weiss 1907) 
Bethe approximation (1935) 


Kirkpatrick and Wolynes 1987 
Sec. 2.1 

Random graphs, cavity method 


Quantitative theory 
in finite d 


Liquid theory (1950s) 
Hypernetted Chain (HNC) 
Percus-Yevick (PY) 


DFT (Stoessel- Wolynes, 1984) 
MCT (Bengtzelius-Gotze-Sjolander 1984) 
Replicas (Mezard-Parisi 1996) 
Jamming - Sec. 2.2 


Corrections 
around MFT 


Ginzburg criterion, d u = 4 (1960) 

Renormalization group (1970s) 
Nucleation theory (Langer, 1960) 


Ginzburg criterion, d u = 8 (2011) 
Renormalization group (2010-) 
Nucleation theory (KTW 1987) 



Table 1.1: A very schematic view of the analogy between RFOT theory and the theory of standard 
phase transitions. References in bold point to later sections where original material will be presented. 



1.2 The random first order transition theory of the glass tran- 
sition 

The random first order transition (RFOT) theory is a very ambitious program to obtain a rather 
complete theory of the liquid-glass transition and of the glass state. It basically consists in following 
the same approach that led to the classical theory of phase transitions, adapting it to take into account 
the presence of disorder. Here a very schematic presentation of RFOT theory will be given, based 
on the analogies of Tab. 1.1. More detailed presentations, and a comparison with other theories of 
glasses, can be found in [41, 14, 18, 42, 43]. 

1.2.1 An overview of the classical theory of phase transitions 
Mean field theory 

The classical theory of phase transitions started from a very schematic theory, the Landau theory [44] . 
Landau theory allows one to identify the correct order parameter, which for the liquid-gas transition 
is a scalar order parameter that signals the onset of spontaneous Z 2 symmetry breaking. It is a 
schematic description that generalizes previous theories such as the Van der Waals equation of state 
of liquids [45] and the Curie- Weiss description of magnetism [46, 47]. These two descriptions have 
also the advantage of being exact in some well defined limit of given microscopic models [48, 49]. In 
particular, they are both exact when the spatial dimension goes to infinity. Take for instance the Ising 
model in d dimensions, with a ferromagnetic coupling J / (2d). For d — > oo, the solution of this model is 
given exactly by the Curie- Weiss equation m = tanh(/3 Jm) [50]. The reason is that when d — > oo, the 
number of neighbors of a given spin grows quickly, and correlations among them are correspondingly 
reduced, because each of them is connected to a large number of other neighbors. At the leading 
order in 1/d, the model on a hypercubic lattice is well approximated by the model defined on a tree 
of the same connectivity, which is known as the Bethe approximation [51] and provides an improved 
version of the Curie- Weiss theory. Similarly, the Van der Waals equation provides an exact theory of 
the liquid-gas transition when d — > oo for liquids with a hard core and a properly scaled attractive 
interaction. Because all these theories are based on the idea that neighbors of a given particle are 
uncorrelated and provide a "mean field" acting on the particle, to be determined sclf-consistcntly, 
they are known as mean field theories (MFT) . 
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Approximate theories in finite dimensions 

The Van der Waals theory can be seen as a simple resummation of the low-density virial expansion. 
Improved theories of the liquid phase (and sometimes also of the liquid-gas transition) can be obtained 
by more careful resummations of the virial expansion. This leads to self-consistent integral equations 
for the correlation function g(r) of the liquid. The most popular of such approximations are the 
Hypernetted Chain (HNC) or Percus-Yevick (PY) closures [6]. Although these theories provide a much 
better description than the mean field theory they retain a mean field flavor as they somehow neglect 
three-body correlations and provide a self-consistent equation for the pair correlations. Improved HNC 
and PY equations, that contain three-body correlations, can be obtained. These theories provide 
reasonable quantitative characterizations of the liquid phase for a given choice of the microscopic 
interaction potential, provided one is not too close to a critical point. 

Corrections around mean field theory 

A somehow different approach focuses on the description of the vicinity of the critical point, where the 
system is characterized by a large correlation length and some of its properties become universal, in 
the sense that they do not depend on the microscopic details but only on the type of symmetry that is 
broken at the transition. It is observed that mean field theory gives good predictions for the universal 
properties for large enough d, but the predictions become poor below some upper critical dimension 
d u . Hence in this case a good approach is to start from the qualitative Landau-Ginzburg theory and 
try to perform a systematic inclusion of non-mean field effects. The simplest calculation that one can 
do is a one loop perturbative calculation. Asking that the one loop corrections remain smaller than the 
leading mean field results leads to the Ginzburg criterion (see [52] and references therein) . The latter 
states that above dimension d u = 4, mean field theory gives the correct critical exponents. Still loop 
corrections will be relevant for non- universal quantities if one is sufficiently close to the critical point. 
How close to the critical point one should be to see deviations from mean field predictions depends on 
the details of the model, but in any case when d — > oo the interval where corrections are important 
shrinks to zero. On the contrary, below d u — 4 loop corrections affect also the universal predictions of 
mean field theory, i.e. the critical exponents. In this case the Ginzburg criterion is however universal: 
it states that deviations from mean field become visible when the correlation length of the system 
is larger than a universal number. In this case the mean field predictions for the critical exponents, 
and in some cases also the predictions for non-universal quantities, can be improved through the 
renormalization group approach [53, 50, 54]. 

Another important direction to improve mean field theory is to add the effect of nuclcation below 
the critical point, where a first order transition is observed. Here, starting from the mean field theory, 
one can take into account the effect of nucleation by non-perturbative instantonic calculations or via 
phenomenological theory [55] . 

Summary 

In summary, the classical approach to phase transitions starts from the qualitative Landau theory, that 
provides a framework to understand the universal features of the transition in terms of spontaneous 
symmetry breaking and the associated order parameter. Then, the theory can be improved in two 
independent directions. On the one hand, one can try to be more quantitative by keeping the mean 
field nature of the theory but taking into account more microscopic details of the system: this is the 
basis of the liquid theory developed in the 1950s, and it works well far enough from the transition. 
On the other hand, one can remain qualitative (i.e. forget about the microscopic details) but try 
to include non-mean field effects: mainly the effect of critical fluctuations at a second order phase 
transitions, that are included perturbatively through the renormalization group, and the effect of 
nuclcation at a first order transition point, that are non-perturbative in nature. This approach is 
presented in Tab. 1.1. 

Of course one would like in the end to develop, for a given system, a complete theory capable 
of describing correctly both its qualitative (universal) properties and quantitative (non-universal) 
properties. However, already for simple liquids this seems to be a formidably difficult task. 
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1.2.2 Qualitative mean field theory of the glass transition 

The RFOT theory aims at repeating the same steps discussed above for standard phase transitions in 
the case of the liquid-glass transition. The first step is then to identify a qualitative, Landau-like theory 
that reproduces the basic facts discussed in Sec. 1.1. Next, one would like to construct quantitative 
approximate theories in finite dimensions, and finally understand non-mean-field corrections related 
to both critical fluctuations and nucleation - as it will turn out that both corrections simultaneously 
contribute in glasses. Before proceeding, it is worth to stress that of course nothing guarantees that the 
same program that is successful for standard phase transitions should work for the glass transition. 
Maybe the latter is a completely different phenomenon that requires an entirely new theoretical 
approach. Therefore, the predictions of the RFOT theory must be, as usual in physics, tested a 
posteriori against experiment and numerical simulations. The result of such tests is encouraging as 
discussed e.g. in [41, 43, 1, 56] and in Sec. 2.2. Of course this does not mean that alternative -and 
often complementary- approaches to the glass transition are less interesting, see e.g. [42, 4] for a 
review. 

The p-spin universality class 

The qualitative mean field theory of the glass transition was constructed by Kirkpatrick, Thirumalai 
and Wolynes (KTW) in a series of pioneering papers appeared from 1987 to 1989 [57, 58, 59, 60, 61, 
62, 63]. Building on the concept of replica symmetry breaking developed by Parisi in the solution of 
the mean field Shcrrington-Kirkpatrick (SK) model of a spin glass [64], KTW realized that a simple 
generalization of the SK model has the desired phenomenology. This is the so-called p-spin glass, 
whose Hamiltonian is given by 



where cr, are either real variables subject to a spherical constraint J2 i of = N, or Ising variables, 
<Tj = ±1, and Ji lt ... t i are independent quenched random Gaussian variables with zero mean and 
variance p\ J 2 /(2N P ~ 1 ). The sum is over all the ordered p-uples of indices i\ < ■ ■ ■ < i p . It is a mean 
field model because each degree of freedom interact with all the others with a strength that vanishes 
in the thermodynamic limit, in order to have an extensive average energy. 

This system shows (Fig. 1.7) an equilibrium Kauzmann transition at a finite temperature Tk, 
where the configurational entropy vanishes, the specific heat jumps downward and the order param- 
eter discontinuously jumps to a finite value. Its dynamics is very similar to the one of glass forming 
liquids in the region of temperature T m > T » T g , but the VFT behavior of the relaxation time 
is not reproduced by these models: instead, a power law divergence of the relaxation time is found 
at a temperature Td > Tk- Although this scaling is due to the mean field nature of this model, 
it is not completely unrelated to what is observed in glass forming liquids, where the behavior of 
T a at temperature T not too close to T g can also be described by a power law. Indeed, the equa- 
tions that describe the dynamics of the p-spin glass model are formally very similar to the so-called 
Mode-Coupling equations [56] that describe well the dynamics of supercooled liquids in a range of 
temperature below T m but not too close to T g [67]. The dynamics displays aging if the system is 
prepared at equilibrium above Td and quenched below this temperature [68] . This model also displays 
dynamical heterogeneities: one can define a dynamical susceptibility of the form of Eq. (1.10), which 
diverges on approaching Td [59, 69]. 

The study of the p-spin model provided a crucial connections between statics and dynamics, 
because the dynamical transition is found to be related to the emergence of metastable states of the 
system below T d , separated by infinite barriers [59, 60, 62]. These metastable states have infinite life 
time in the mean field framework, and they are responsible for the divergence of the relaxation time 
at Td- It is important at this point to anticipate that in finite dimension, it is reasonable to expect 
that the life time of the metastable states will be finite. Activated processes of barrier-crossing will 
restore the crgodicity below Td, and the relaxation time will only diverge at Tk- A relatively simple 
nucleation theory shows that the relaxation time should diverge in a VFT-like way at T K [61], and 
predicts the existence of a static divergent correlation [23, 25]. 
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Figure 1.7: Qualitative mean field phase diagram of RFOT theory, as described in Sec. 1.2.2. A 
dynamical transition is observed at Td, where r Q ~ (T — X^) -7 . Below Td, the configurational entropy 
is finite and many glasses dominate the thermodynamics, while the dynamics is not ergodic. At Tk, 
the configurational entropy vanishes and the number of glasses becomes finite. It is reasonable to 
conjecture that in finite dimensional systems activated processes restore ergodicity below Td so that 
Tq, only diverges at Tr- according to a VFT-like law (red dashed line). The same phase diagram can 
be obtained as a function of increasing density instead of decreasing temperature in Bethe lattice glass 
models [65, 66]. 

Although p-spin models seem quite abstract and very different from structural glasses, one can 
construct mean field models that belong to the same universality class and are formulated in terms of 
interacting particles [65, 66]. This is interesting because these particle models also display a. jamming 
transition with a phase diagram identical to the one in Fig. 1.6. Moreover, these models do not display 
any quenched disorder, like structural glasses, hence suggesting that the presence of quenched disorder 
in the p-spin Hamiltonian is not an essential ingredient to produce a glassy phenomenology. 

The Landau theory that describes the glass transition of the p-spin universality class is formulated, 
by using the replica trick, in terms of an overlap matrix q a b, which is the order parameter of the 
transition. It is the self-overlap between two different identical copies (replicas) of the system, labeled 
a and b: 



and plays the role, in the context of spin glass theory, of the nonergodicity factor (1.5). If q a b(x) is a 
local version of the overlap in point x, the Landau free energy has the form [70, 71, 72, 73] 



For the moment the number of replicas is left unspecified, this point will be discussed below. 

In summary, the basic facts listed in Sec. 1.1.6 are reproduced by this mean field theory. Moreover, 
the theory is predictive: the existence of dynamical heterogeneities [59, 69] and of a static correlation 
length [23, 25] were first predicted based on the analysis of p-spin-like models, and later observed 
numerically. Excellent reviews on the properties of the p-spin-like models have been recently pub- 
lished [74, 75, 67]; in the following only the main results will be reviewed, referring to these works for 
all the details. 
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The TAP free energy 

To better understand what is going on in these models one has to investigate the structure of their 
phase space. In particular, one wishes to characterize the equilibrium states in order to understand 
the nature of the thermodynamical transition at Tk, as well as the structure of the metastable states 
that trap the system at Td and are responsible for the existence of a dynamical transition. 

It turns out that at the mean field level a pure state a (equilibrium or metastable) is completely 
determined by the set of local magnetizations mf, i = 1, • • • ,N. The local magnetizations are minima 
of the Thouless- Anderson-Palmer (TAP) free energy [76, 77] which is defined by 



F(rm) = -riog I Y e- pH W +f3 E* M^-™*) L (1.15) 



where the auxiliary fields hi are introduced to fix the local magnetizations and are determined by the 
condition dF/dhi = at fixed . Its explicit calculation relies on high temperature or large dimension 
expansions [78, 79]. The weight w a of state a is proportional to exp[— f3Nf a ], where f a — F(mf)/N. 
Local minima of F having a free energy density / > f m in f° r N — > 00 are metastable states. The TAP 
free energy F{rrii) depends, in general, explicitly on the temperature, so the whole structure of the 
states may depend strongly on temperature. 

In mean field p-spin models, at low enough temperature, the number of states of given free energy 
density / is £!(/) = exp[iVX(/)]. The function £(/) is a decreasing function of /, that vanishes 
continuously at / = f m in and jumps discontinuously to zero above / = f max - A similar behavior is 
found in all models belonging to the p-spin universality class. The main peculiarity of p-spin models 
is that an exponential number of metastable states is present at low enough temperature. 

One can write the partition function Z, at low enough temperature and for N — > 00, in the 
following way: 

z = e -fiNF{T) „ J- e -0N fa = [ }maX df e N[HU)-m „ e N[Z(n-0fl , (L16) 

Jfmin 

where /* <G [fmin, fmax] is such that / — TT,(f) is minimum, i.e. it is the solution of 



effi _ 1 
W ~ T 



(1.17) 



provided that it belongs to the interval [f m in, fmax]- Starting from high temperature, one encounters 
three temperature regions (Fig. 1.7): 

• For T > Td, the free energy density of the paramagnetic state is smaller than / — TS(/) for any 
/ £ [fmin, fmax], so the paramagnetic state dominates and coincides with the Gibbs state (in 
this region the decomposition (1.16) is meaningless). 

• For T d >T > T K , a value /* g [f m in, fmax] is found, such that /* - TS(/*) is equal to f para . 
This means that the paramagnetic state is obtained from the superposition of an exponential 
number of pure states of higher individual free energy density /* . The Gibbs measure is split on 
this exponential number of contributions: however, no phase transition happens at Td because of 
the equality /* — TS(/*) = f para which guarantees that the free energy is analytic on crossing T^. 

• For T < Tk, the partition function is dominated by the lowest free energy states, /* = f m in, 
with S(/ min ) = and F(T) = f min - TT,(f min ) = f min . At T K a phase transition occurs, cor- 
responding to the 1-step replica symmetry breaking transition found in the replica computation. 
The order parameter q a b jumps discontinuously at the transition. 

In the range of temperatures T^ > T > Tk, the phase space of the model is disconnected in an 
exponentially large number of states, giving a contribution S(T) = E(/*(T)) to the total entropy of 
the system. This means that the entropy S(T) for Td > T > Tk can be written as 

S(T) = E(T) + S mb (T) , (1-18) 
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SvibiT) being the individual entropy of a state of free energy /*. From the latter relation it turns out 
that the complexity S(T) is the p-spin analogue of the configurational entropy S C (T) of supercooled 
liquids 1 . 

The TAP approach provides also a pictorial explanation of the presence of a dynamical transition 
at Td- If the system is equilibrated at high temperature in the paramagnetic phase, and suddenly 
quenched below Td, the energy density start to decrease toward its equilibrium value. This relaxation 
process can be represented as a descent in the free energy landscape at fixed temperature starting 
from high values of /. When the system reaches the value f max it becomes trapped in the highest 
metastable state and is unable to relax to the equilibrium states of free energy /* , as the free energy 
barriers between different states cannot be crossed in mean field [75, 67]. For this reason below Td 
the system is unable to equilibrate. 



Real replica method 

For systems that belong to the p-spin universality class, two very powerful and general methods to 
compute the complexity as a function of the free energy of the states without directly solving the 
TAP equations exist; they have been developed in [60, 87, 88, 89, 90, 91, 92, 93]. They go under the 
name of the real replica method [89] and the potential method [88]. Both methods consider a number 
of copies of the system coupled by a small field conjugated to the order parameter q^. Here, only a 
very short account of the real replica method will be presented. 

The idea of [60, 89, 87] is to consider m copies of the original system, coupled by a small attractive 
term added to the Hamiltonian. The coupling is then switched off after the thermodynamic limit has 
been taken. For T < T d , the small attractive coupling is enough to constrain the m copies to be in 
the same TAP state. At low temperatures, the partition function of the replicated system is then 

/•/ma* 

0Nt>(m,T) _ ^ e -[3Nmf a = l d j e JV[E(/)-/3m/] _ p JV[S(/* )-/3m/*] 

a J fmin 



Z m = e - p ^ m ' I >~ Y, e ~ P ta = / dfe N W>- pm " ~ } - pm} 1 , (1.19) 

where now f*(m,T) is such that mf — TE(/) is minimum and satisfies the equation ^ = ?p. If m 
is allowed to assume real values by an analytical continuation, the complexity can be computed from 
the knowledge of the function $(m,T) = mf*(m,T) — TS(/*(m,T)). Indeed, it is easy to show that 

a<i>(m,r) 

/ ( m ' T ) = — — . 

1 (1.20) 
S(m, T) = £(/*(m,T)) = m 2 ^ m ^ m,T ^ = m/3f*(m,T) - 0$(ro,T) . 

The function £(/) can be reconstructed from the parametric plot of f*(m, T) and E(m, T) by varying 
m at fixed temperature. In other words, $(m, T) is the Legendre transform of the function £(/). 
Its knowledge allows one to obtain the function £(/) and from it reconstruct all the thermodynamic 
properties of the system, including those of metastable states and the dynamical transition Td- In [92] 
this method was applied to the spherical p-spin system and it was shown that the method reproduces 
all the results obtained from the explicit TAP computation. Note that replicas are introduced here 
to take into account the existence of multiple equilibrium states, without any reference to quenched 
disorder. 



Dynamics 

The Langevin dynamics of the model can also be solved exactly [59, 68, 67]. The Martin-Siggia-Rose 
formalism [94, 95] allows one to write it in terms of a closed equation for the time-dependent spin-spin 

1 h\ the interpretation of experimental data one should remember that in experiments S v n, can be estimated only 
by the entropy of the crystal. However, the vibrational properties of the crystal can be different from the vibrational 
properties of an amorphous glass, sec [80] for a review. Corrections due to this fact must be taken into account: in some 
cases it has been claimed that the difference is reduced to a proportionality factor between S c and E [81, 82]. On the 
contrary, in numerical simulations £(T) and S v n,(T) can be measured directly [83, 84, 85, 86]. 
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correlation function 

iV— >oo 1\ ^— ' 

Note that the latter measures the overlap between the configuration at time and that at time t. The 
self-consistent equation for C{t) reads [67, 74] 



eft 



= -TC(t)-£f od uC(t-uy-^. (1.22) 



From this equation, one can show that a dynamical transition is found at the same temperature T d 
that comes from the TAP computation. Upon approaching T d , C(t) has the same form as in Fig. 1.2. 
It shows a two-step relaxation, with well separated time scales and a structural relaxation time that 
shows a power-law divergence for T — > T d , r a <~ |T — T d |~ 7 . A dynamical order parameter can be 
defined as 

% = lim C(t) ; (1.23) 

t— >oo 

it is the analogue of the dynamical nonergodicity factor defined in (1.5) and jumps to a finite value 
at T d . Below T d the system is no more able to equilibrate with the thermal bath and enters a 
nonequilibrium regime. This result confirms that metastable states are responsible for the slowing 
down of the dynamics and for the dynamical transition at T d . 

As already mentioned, a four-point dynamical susceptibility associated to C(t) can be defined, and 
one can show that it diverges at T d [59, 69]. 



1.2.3 Quantitative mean field theory (theories) of the glass transition 

Even before the coherent RFOT theory was formulated by KTW, several independent attempts to 
construct quantitative approximate theories, in the same spirit of the HNC or PY theory of liquids, 
had been formulated. Here a short review of these theories will be given, focusing in particular on 
their connection with the broader RFOT picture. 



Density functional theory 

A density functional theory (DFT) of glassy states was first proposed in [96, 97] and later improved 
by several authors [59, 60, 98, 99, 100, 101]. DFT is the analog of the TAP approach for particle 
systems. One writes the free energy of the system as a functional of the density field. A commonly 
used approximate form is 

-pF[p(x)] = J p(x)[l - logp(x)} + l -J dxdx'p(x)p(x')c(x - x 1 ) , (1.24) 

where c(x) is the direct correlation function of the liquid [6]. The density profile is determined 
by minimization of the free energy. Because in general this is a too difficult task, one looks for a 
parametrization of the density field as a sum of Gaussian terms: 

AT (^--"j) 2 

P{x) = £ e —^r ■ (1.25) 
j=i ^2^,4 

Now the parameters A (the "cage radius") and Rj (the "equilibrium positions" in the glass) have 
to be determined by inserting this parametrization of the density in Eq. (1.24) and minimizing the 
result, which has the form 

-PF[p(x))/N = d - \og{2-KA) + \ J ^c^So^e-^ 2 , (1.26) 
where Sa(q) = N^ 1 J2j k e lq ^ Rj ~ Rk ^ is the structure factor of the equilibrium positions. 
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At this point one can either choose to approximate So(q) with the liquid structure factor, So(q) ~ 
S(q), and determine only A by minimization [96, 97] or try to perform a full minimization including 
the Rj for a finite system size N [98, 99, 100, 101]. In both cases, amorphous solutions with a finite 
A are found at high enough density or low enough temperature. However, DFT has several problems: 

• The approximations that are made, starting from the definition of the density functional in 
Eq. (1.24), are quite crude and for this reason the quantitative accuracy of DFT is very poor. 
Of course one could improve these approximations, but a satisfactory quantitative DFT has not 
been obtained for the moment. 

• Even if some amorphous solutions can be found numerically at finite N, counting their number 
and its scaling with A^ is a formidable task, hence a proper thermodynamic characterization of 
the glassy states and their configurational entropy £(/) is impossible. 

For these reasons, DFT remains for the moment only useful to illustrate the basic mechanism of 
the glass transition in the RFOT framework: the discontinuous appearance of amorphous free energy 
minima [96, 97, 59, 60]. A satisfactory quantitative description of this phenomenon in d — 3 requires 
considerable improvement of DFT. 

Mode-coupling theory 

Mode-coupling theory (MCT) has been initially developed to describe the critical slowing down at a 
second order phase transition [102], and later generalized to describe the glass transition [103, 56]. 
MCT equations are closed equations for the coherent scattering function defined in Eq. (1.4). If 
<f>(k,t) = F(k,t)/S(k) is the normalized correlator such that <f>(k,t = 0) = 1, the MCT equations 
read [104, 56] 



where fid is the (i-dimensional solid angle and 7^ is a /c-dependent microscopic relaxation rate related 
to S(k). The kernel V(k,p, q) is a function of the static liquid correlation functions S(q) and c(q) [104]. 

The analogy between MCT equations (1.27) and the exact equations that describe the Langevin 
dynamics of p-spin models, Eq. (1.22), is striking. In particular choosing p = 3 one sees that the only 
difference is the fact that the MCT equations are complicated by the g-dependence of all correlation 
functions, that is clearly due to the finite dimensional nature of MCT. This similarity has led many 
to conjecture that MCT is the finite-dimensional formulation of the dynamic theory of RFOT [62, 59, 
60, 75, 67]. However, for the moment, the MCT equations have always been derived by performing 
some "uncontrolled" approximations, in the sense that no systematic improvement of the derivation 
can be obtained. A derivation of MCT in terms of a systematic expansion (possibly a low-density 
or high-dimensional expansion) [105, 106, 107] would help clarify its relation with the general RFOT 
scenario. It has been shown that the basic phenomenology of MCT is insensitive to the detailed form 
of the kernel, suggesting that indeed MCT can be seen as a Landau theory of the dynamics at a glass 
transition [108]. 

Despite the fact that the direct connection between MCT and the RFOT picture remains partially 
unclear, the MCT equations can be solved both numerically and analytically in some asymptotic 
regimes and their predictions are identical to those obtained from the mean field Eq. (1.22): they 
predict a sharp discontinuous transition at a temperature T<j, where the relaxation time of 4>(k,t) 
diverges as a power-law. Below T<j, the non-ergodic parameter lim t _ +00 <j>(k, t) = f(k) and the cage 
radius A are both non-zero as in Eq. (1.5). MCT gives quantitative predictions for Td and f(k). 
Moreover on approaching Td from above, it predicts a number of scaling laws. Detailed tests of these 
predictions have been made both numerically [10, 11] and experimentally [109, 110, 111]. There is 
quite general consensus that MCT predictions work very well in a range of temperatures slightly 




(1.27) 




k+p 



18 



CHAPTER 1. A THEORY OF THE GLASS TRANSITION 



above the glass transition T gi while on approaching T g deviations from MCT behavior become more 
pronounced. 

An important point is that MCT did not limit itself to provide quantitative results for already 
known quantities. It also provides new predictions for previously unknown phenomena, that were 
later confirmed by numerical simulations and experiments. The most impressive one is probably the 
prediction of a re-entrant glass transition line and the existence of a glass-liquid-glass and a direct 
glass-glass transition in attractive colloids [112]. 



Replica theory 



Despite early attempts [59] , a quantitative replica theory (RT) of finite dimensional glasses has been 
developed much later than DFT and MCT, mainly thanks to the work of Mezard and Parisi [113, 114] 
who found a suitable formulation of the real replica method of Monasson [89] for particle systems. 

The dynamical transition can be described in a static framework by introducing a replicated version 
of the system [89, 113]: for every particle one introduces m — 1 additional particles identical to the 
first one. In this way one obtains m copies of the original system, labeled by a — 1, . . . , m. The 
interaction potential between two particles belonging to replicas a, b is v ab (r), where v aa (r) = v(r), 
the original potential, and v ab (r) for a ^ b is an attractive potential that constrains the replicas to be 
in the same metastable state. The basic fields are the one and two point density functions 



N 



Pa{x) =J2 S ( X ~ X i) 



(1.28) 



P { ab'( X >y) = Pa{x)Pb{y) ~ pa{x)S ab 5(x - y) 



To detect the dynamical transition one has to study the two point correlation functions when 
v a b{f) — > for a 7^ 6, and in the limit m — > 1 which reproduces the original model [89, 113]. In this 
limit, the two-replica correlation function is, for a ^ b: 



(1.29) 



(C ab (r)) = / dxf(x) (p a (r + |) p b (r - |) ) - p 2 . 



Because of the limit v ab (r) — > 0, the two replicas fall in the same state but are otherwise uncorrelated 
inside the state, therefore (C ab (r)) = (C(r,t — > oo)) which provides the crucial identification between 
replicas and dynamics [73, 1]. Similar mappings can be obtained for four-point correlations. 

One can introduce for convenience an external field v a (x) (that is derived from a space-dependent 
chemical potential), in such a way that the density correlation functions can be obtained by tak- 
ing the derivative of the free-energy with respect to it [6]. The free energy is defined as the log- 
arithm of the partition function, and its double Legendre transform defines the Gibbs free energy 
n{Pa(x)},{p$(x,v)}] [6, 115]: 



dxdy 



Pab( x >y) ln 



Pab( X >V) 



- Pab( x >y) + pa{x)p b {y) 
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n>3,ai,...,a„ 



2n 



Pa{x)p b {y) 

dxi ■ ■ -dx n p ai (xi)h aia2 (xi,X 2 ) ■ ■ ■ Pa n ( x n)ha nai ( x n,Xl) + T 2 PI 



+ E / dx P^ X ) \ Xn Pa{x) - 1] 



(1.30) 



(2) 

where h ab (x,y) — p ab (x,y)/ p a {x)p b {y) — 1 and T2Pi is the sum of 2-line irreducible diagrams [115]. 
The average values of the fields in Eq. (1.28), namely p a (x) and p ab (x,y), can be obtained by solving 
the saddle point equations 



STiiPajAP^}} 



s Pab( x >y) 



= ^v ab (x,y) 



(1.31) 



and similarly for p a {x). 
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The crucial advantage of replica theory with respect to DFT is that the whole system of m replicas 
is in its liquid phase, hence it is homogeneous. In fact, introducing replicas allows one to sum over all 
amorphous metastable glasses, thus recovering translational invariance. Therefore, one can restrict 
the considered solutions to those where p a {x) = p, a choice that allows one to greatly simplify the 
computation with respect to DFT, where instead one has to work with the full amorphous density 
profile p(x). 

A concrete implementation of this procedure amounts to neglecting the 2PI diagrams, therefore 
obtaining a replicated version of the HNC equations [113]. In this way one obtains predictions for T4, 
Tk, the complexity function, and the non-ergodic parameter f(k) [113]. The results for Td and Tk 
are reasonably good, however the prediction for f(k) is much worse than the MCT one. Moreover, 
the complexity is found to be quite a bit smaller than the numerical estimate, and this scheme gives 
inconsistent predictions in the glass phase below Tk [113]. 

A better approximation scheme in replica theory is obtained by considering a slightly different 
formulation, in which the basic degrees of freedom are molecules built of an atom from each replica. 
The physical reason for this is that replicas are assumed to be in the same state, hence the density 
profile of the different replicas is assumed to be similar, and this can be described by assuming that 
atom 1 in replica 1 is very close to atom 1 in replica 2, and so on. Furthermore, the average distance 
between atoms in the molecule, which is related to the cage radius A, is assumed to be small and one 
performs a systematic expansion in powers of A starting from the standard formulation of molecular 
liquid theory (e.g. in the HNC approximation) [87, 114]. This formulation is much more effective for 
treating the glass phase at Tk and below, where the cages are well formed. Indeed, this approach 
gives good quantitative predictions for the configurational entropy above Tk, and for the specific heat 
and the energy of the glass phase below Tk [87, 84]. Moreover, this approach allows one to obtain a 
very detailed theory of the jamming transition that happens deep in the glass phase. This will be the 
subject of Sec. 2.2. The main drawback of this approach is that because molecules are assumed to be 
stable and have small radius, the process of "molecule dissociation" that in replica language leads to 
the disappearance of metastable states at Td is not very well described by this approach. Depending 
on the details, either Td is not found, or it is largely overestimated. 

The construction of a "mixed" replica approach, capable of describing both the region around and 
below Tk where cages are well formed, and the region around Td where cages are unstable, would 
be an important achievement, that could be obtained by mixing the diagrams of the replicated HNC 
with those of the molecular liquid formulation. 

Towards a consistent finite dimensional implementation of RFOT theory 

The different aspects of the RFOT scenario for mean-field p-spin glasses all find a counterpart in the 
study of the liquid-glass transition in finite dimension. The equivalent of the TAP approach is DFT: in 
both cases one tries to find the metastable states by minimizing a suitable free energy functional. The 
dynamics of p-spin models is described by an equation that corresponds to the so-called "schematic" 
limit of MCT, i.e. the limit where the wavevector dependence is neglected. Finally, replica theory can 
be formulated for p-spin models and liquids along the same lines. 

Yet, the important difference is that in mean field p-spin models, it can be shown exactly that the 
TAP, dynamic and replica approaches give completely equivalent descriptions of the same phenomenon, 
which is the emergence of glassy metastable states [74]. While in finite dimensions, DFT, MCT and 
RT are necessarily approximate theories that use different approximations and lead to quite different 
quantitative predictions, with marked discrepancies e.g. on the value of Td, on the shape of the 
non-ergodic parameter, and so on. This is not particularly worrying as far as d = 3 is concerned, 
because already in liquid theory different approximations schemes (e.g. HNC and PY) obviously lead 
to different predictions [6]. 

However, for the RFOT scenario to be consistent, one would like to show that in the limit d — > 00, 
the static and dynamic descriptions of the transition become equivalent, exactly in the same way as the 
HNC and PY equations both converge to the Van der Waals equation of state. This consistency check 
was attempted in one of the earliest works on RFOT theory [62]. There it was shown that under a 
Gaussian approximation for the non-ergodic parameter, DFT and MCT lead to equivalent predictions 
for the dynamical transition point of hard spheres. Another interesting attempt in this direction is the 
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work of Szamel who derived the MCT equations from a static replica approach [116]. However, despite 
these apparent successes, the limit d — > oo turns out to be much more complicated: the solution of 
full MCT (without the Gaussian approximation) [104, 117] leads to a strong discrepancy with RT [1] 
and DFT [62]. Therefore despite these results, the problem of obtaining a consistent RFOT theory of 
the glass transition in d — > oo remains open. This will be discussed in Sec. 2.1. 

1.2.4 Corrections around mean field theory 

After a consistent mean field theory of the glass transition has been obtained, in finite dimensions one 
would like to compute the corrections due to fluctuations around the mean field approximation. When 
this program is carried out, one finds that (according to the mean field theory itself) there are two 
important sources of corrections to the mean field scenario. The first corrections originate from critical 
fluctuations (the dynamical heterogeneities) that become important around the dynamical transition 
Td below the upper critical dimension, as in any standard critical phenomenon [118, 73]. The second 
corrections are non-perturbative phenomena related to activated processes, that are important below 
Td, close to the experimental glass transition T g and down to Tk- They can be taken into account 
by a phenomcnological approach, leading to a number of predictions that are in good agreement with 
experiment [41]. This approach can be partially justified via a renormalization group analysis [119, 72]: 
however, this part of the theory remains the most controversial [120]. 

One big source of controversy is the very existence of the Kauzmann temperature Tk, which relies 
on a dubious extrapolation of the relaxation time below T g . In fact, many other functional forms 
that do not display a singularity at finite temperature are compatible with the data, see e.g. [17]. 
Many theoretical arguments have indeed been proposed against its existence in finite dimensions [121, 
122]. None of them is however completely convincing, leaving the question open to debate. The 
disappearance of the Kauzmann transition in any finite dimensions (or possibly below some lower 
critical dimension larger than 3) could be due to a non-perturbative process that is not taken into 
account in the mean field theory. For this reason, the present approach (start from MFT and then 
consider corrections around it) seems unable to capture such hypothetical processes. In any case, 
although the existence or not of Tk is certainly an important question from the theoretical point of 
view, it is not so important for the theory of everyday glasses. In these systems, one is always quite far 
from Tk, and the static and dynamic length scales never reach too large values (actually, they typically 
are of the order of a few interparticle distances at T g ). Therefore, the really important question for 
the theory of glasses is not whether MFT retains its validity at all length scales, down to the transition 
point. The important question is whether MFT (plus corrections) retains its validity down to T g , and 
whether it can explain the phenomenology of the glasses we observe in nature. It might well be that 
MFT provides a good description in a region where the relevant length scales are smaller than, say, 
100 interparticle distances, and for larger length scales some unknown process completely washes away 
the mean field phenomenology. If this is the case, MFT would be an excellent theory of real glasses, 
because one never observes such large length scales and deviations from mean field would remain 
unobserved. Such a situation would be somehow similar to the BCS theory of superconductivity, 
which is a mean field theory but yet provides excellent qualitative and quantitative predictions. This 
is due to the fact that the pairing length a is so large that each Cooper pair effectively interact with 
a large number of other pairs. Hence the Cooper pair correlation length £ is much smaller than a, 
£/a <C 1, making MFT a very good description. Indeed, the associated Ginzburg criterion tells that 
one has to be extremely close to the transition to observe deviations from BCS theory. The aim of 
the rest of this discussion will be therefore to give criteria of validity of MFT sclf-consistently within 
MFT itself. 

Critical fluctuations around MCT 

Within the static formulation of RFOT theory, the dynamical transition at T<j is a spinodal point at 
which the metastable states become unstable and disappear [73, 123]. Associated with the instability 
one finds a diverging correlation length (akin to the diverging length that is found at a spinodal point) 
that gives rise to dynamical heterogeneity [69]. As for a standard spinodal point, in its vicinity one 
can reformulate the field theory in Eq. (1.14) as a cubic theory. It has been shown in [73] that the 
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relevant theory is a ip 3 theory in a random field. Performing a systematic loop expansion around 
the mean field saddle point in such a theory, one can obtain a Ginzburg criterion and compute the 
upper critical dimension d u . This allows one to establish that d u = 8 (for the (p 3 theory d u = 6, but 
the random field provides dimensional reduction) [73]. Using the HNC replicated free energy given 
in Eq. (1.30), a quantitative estimate of the Ginzburg criterion has been given in [124] for d < d u . 
The result of this analysis is that in d = 3, one needs a dynamical correlation length of order 1 (in 
units of interparticle distance) to observe deviation from the mean field behavior. A comparison with 
numerical data obtained in [125] shows that the dynamical correlation becomes of order 1 only slightly 
below Td- Hence in the full range where the MCT phenomenology is observed, the Ginzburg criterion 
indicates that non-mean-field corrections remain small. Below Td, these corrections might become 
important, however at that point activated processes start to dominate over the MCT phenomenology. 
One concludes that in d = 3 critical fluctuations around the dynamical transition do not play any 
observable role. 

Such critical fluctuations have been also described within MCT [126, 118, 127, 128], leading to a 
very detailed characterization of the dynamical heterogeneity and to similar results for what concerns 
the range of validity of mean field theory In particular, in [120] a phcnomenological Ginzburg criterion 
for MCT has been proposed, that leads to the same conclusion as the analysis of [124]. 

These results seem to indicate that developing a systematic RG theory to discuss critical fluctu- 
ations around the dynamic transition might not be worth the effort, because deviations from mean 
field should be very difficult to observe at least in d = 3. The situation might be different in higher 
dimensions where nucleation is suppressed and the MCT regime is more clearly observed [129]. 

Activated processes and nucleation theory at low temperatures 

The last, and actually most important, ingredient of the RFOT scenario is a phenomenological treat- 
ment of the activated processes that restore ergodicity between Td and Tk- These processes are at the 
origin of the slowing down at T g and in particular they are responsible for the exponential divergence 
of the relaxation time, Eq. (1.1). In the MFT, the transition at Tk has a mixed nature. It is a first 
order transition from the point of view of the order parameter: in fact, at Tk the free energy of two 
phases -the one where replicas are uncorrelated and the one where replicas are in the same state- 
cross. However, it is a second order transition from the thermodynamic point of view. This mixed 
character makes the Kauzmann transition a very peculiar one. Indeed, its description requires setting 
up a nucleation theory [61, 41, 43, 23, 120] -due to the first order nature of the mean field transition- 
but at the same time it predicts the existence of a diverging correlation length as in second order 
transitions, with associated critical exponents. In the following a short account of this theory will be 
presented. 

As usual, in mean field theory nucleation processes -where a nucleus of a phase forms inside 
another phase, paying a surface tension cost but gaining in bulk free energy- are not taken into 
account, because in infinite dimensions the surface and volume terms have the same scaling (L d is the 
same as L d ~ x when d — > oo). Hence, in mean field theory metastable states with an intensive free 
energy higher than the free energy of the ground states, / > f m in, still have an infinite life time. These 
states are responsible for the existence of a finite complexity. Their lifetime is infinite, so they are 
able to trap the system below T d . This is the reason why the dynamical transition, i.e. the divergence 
of the structural relaxation time, happens at a temperature Td > Tk- On the contrary, in a model 
with short range interactions, metastable states have a finite lifetime due to the nucleation of bubbles 
of the stable phase inside the metastable one, so they are not thermodynamically stable. One should 
expect the existence of well defined states with / > / TO , n to be impossible; but the analogy between 
mean field models and real glasses is based on the analogy between the complexity S(T) and the 
configurational entropy S C (T). How can one explain the existence of a finite configurational entropy, 
related to well defined metastable states, in a short range system? 

Moreover, the observed crossover of the relaxation time from a power-law behavior to a Vogel- 
Fulcher-Tamman law (1.1), that happens around Td, is not explained by the mean field theory, which 
predicts a strict power-law divergence of r a for T — > T^. The observation of a finite relaxation 
time below T d is again related to the finite lifetime of metastable states. The system, instead of being 
trapped forever into a state, is able to escape, due to nucleation processes; it is then trapped by another 
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state, and so on. These processes of jump between metastable states are activated processes: the 
system has to cross some free energy barrier in order to jump from one state to another. The relaxation 
time is then expected to scale as T a (T) ~ T^exp \(3AF(T)] , AF(T) being the typical free energy 
barrier that the system has to cross at temperature T. The VFT law and the observation that T <~ Tk 
suggest that the barrier should diverge at T K , AF(T) ~ (T — Tk) -1 ; more generally, the Adam-Gibbs 
formula [13] relates this divergence to the vanishing configurational entropy, AF(T) <~ iS c (T) _1 . It 
is then essential to understand what is really the meaning of S C (T) in finite dimension and why it is 
related to the free energy barrier for nucleation. 

A crucial observation (actually, a theorem [24]) is that the divergence of the relaxation time at 
Tk, in short range systems, is possible only if the cooperative processes of structural rearrangement 
involve atoms that are correlated on a typical length scale £, which diverges at T K . If no divergent 
length scale is present in the system, it is always possible to divide it in finite subsystems, each one 
relaxing independently of the others: and the relaxation of a finite system is bound to happen in 
finite time, if the interactions are finite and have short range. A simple idea that follows from the 
above observation and can explain how the mean field picture is modified in short range systems is 
the following [61, 130, 41, 23, 70, 71]. There exists a typical length scale £(T) over which structural 
relaxation processes happen. If one looks at smaller length scales, the system behaves as if it were 
mean field: metastable states are stable for I < £, yielding a finite local complexity. However, for 
large scales I > £, metastability is destroyed and only the lowest free energy states are stable. For 
T — > T^, £ — » oo, so below Tk a stable ideal glass phase is possible. This idea leads naturally to the 
identification of the configurational entropy S C (T) with the local complexity S(T), and to a derivation 
of an Adam-Gibbs-like relation between the relaxation time r a and £(T). 

These considerations have been formalized in [23] where a concrete procedure to define and mea- 
sure £ has been proposed, and later generalized to improved geometries [30]. Following this work, 
numerical measurements of £ have been performed and have revealed a good agreement with the 
theoretical predictions [27, 28, 131, 29]. Analytic calculations in simplified models with Kac-like long 
range interactions [70, 71, 25] and renormalization group calculations [119, 72] have been performed 
confirming this scenario. 

1.2.5 Summary 

To summarize, the RFOT theory of the glass transition is constructed by following closely the path 
that led to the standard theory of critical phenomena, as sketched in Tab. 1.1. The theory was started 
by the pioneering works of Parisi in 1979 [64], who developed the concept of spontaneous replica 
symmetry breaking and identified its consequences in terms of proliferation of metastable states [76] , 
and of Kirkpatrick-Thirumalai-Wolynes [57, 58, 59, 60, 61, 62, 63], who identified a concrete class of 
simple mean field models that reproduced many basic facts of the phenomenology of glasses. At the 
same time, early attempts to construct approximate quantitative theories of the glass transition, such 
as density functional theory [96], mode-coupling theory [103], and replica theory [113], fit well under 
the RFOT umbrella, and can be thought as different implementations of RFOT. Early consistency 
checks [62] corroborated this interpretation, even if many details of it remain open and a complete 
RFOT theory of infinite dimensional liquid is still lacking (Sec. 2.1). Finally, RFOT provides the 
framework to construct a mean field theory of the jamming transition that happens inside the glass 
phase (Sec. 1.1.5) [65, 132, 66, 1, 2]. This will be the subject of Sec. 2.2. 

Nowadays, most of the research on RFOT theory is focused on a detailed understanding of the 
corrections around mean field theory. There are at least two different sources of corrections. Around 
T^, corrections are due to standard critical fluctuations and can be treated via loop expansions both 
in the statics and in the dynamics, leading to a Ginzburg criterion that states that the upper critical 
dimension is d u — 8, and that in d = 3 non-mean-field effects (of this kind) should be barely observ- 
able [118, 120, 73, 124]. Mean field theory predicts in this region the existence of large dynamical 
heterogeneities and an increasing dynamical correlation length, that are indeed observed [18]. The 
associated critical exponents should be the mean field ones for d > 8 [126, 118, 127, 128] and close to 
the mean-field ones for d < 8, but this for the moment remains controversial [125, 133, 134]. Below T& 
and in low dimensions, activation is the main source of corrections around mean field, and a nucleation 
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theory has to be developed to take it into account [61]. Doing this results in a number of predictions 
about the existence of a static correlation length [130, 23, 70, 71, 24, 41, 30] that are being extensively 
tested in simulation [27, 28, 131, 29]. 

A crucial issue is of course how much of the mean field picture survives in physical dimensions 
d = 2,3. It is very likely that the glass transition does not exist in two dimensions [135, 122], 
suggesting a lower critical dimension di > 2. Whether di < 3 or not, and what is the important 
physical mechanism that determines di, remain a largely open problem. Preliminary renormalization 
group investigations [72, 119, 136, 137] give contradictory results. The crucial difficulty, which has 
been extensively discussed in [138], is the following. There are many simple lattice models that 
display a RFOT at the fully-connected mean field level. However, for most of those, the transition 
either disappears or becomes a continuous spin glass transition as soon as some spatial fluctuations 
are included, e.g. by considering the model on a Bcthe lattice. This is a well known effect: the 
inclusion of fluctuations lowers the barriers between states (for example, it lowers considerably the 
critical temperature in the ferromagnetic Ising model), in such a way that models that display the 
transition at the mean field level can become paramagnetic once fluctuations are included. Of course, 
there are also lattice models that have a robust RFOT, such as the lattice glass model of [65] , however 
they are more complicated to study using renormalization group methods. Hence, the problem of 
finding a simple enough model for which the RFOT is robust enough has not been solved for the 
moment [138]. 

It is important to stress, however, that even if the lower critical dimension is bigger than 3, still the 
non-mean-field physics might only be visible close enough to the critical point (or more precisely, when 
the correlation length is large enough). In that case, mean field theory would still provide a rather 
good description of experimental and numerical systems, whose correlation lengths are never very 
large. The fact that many predictions of mean field theory compare well with numerical simulations, 
as will be shown in Chap. 2, seems to support this point of view. 
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Chapter 2 



Glass and jamming transitions of 
hard spheres 

Here the application of the RFOT theory to the glass and jamming transition of hard spheres will be 
presented. This line of research started in 2005 with the paper [132], where the small cage expansion 
within replica theory previously developed in [87] was adapted to hard spheres. Technically, this was 
a non-trivial task because the small cage expansion of [87] is based on a harmonic expansion and 
therefore fails badly for singular potentials. After this extension was achieved, the work evolved in 
two main directions. The first is the study of the glass transition of hard spheres in dimension larger 
than 3, motivated by the desire to achieve a better understanding of how the mean field limit d — >■ oo is 
approached. The second is the study of the jamming transition that happens deep in the glass phase. 
As discussed in Sec. 1.1.5, the jamming transition is characterized by marked structural changes, 
that arc clearly detected in simulation and experiments, and therefore -besides being an interesting 
phenomenon in itself- it provides a unique opportunity to test any theory of the glass state. 

All the technical details of this work are contained in [1, 2]. Here only the main results will be 
discussed. Note that, in the following, v(x) is the harmonic sphere potential, using the same notations 
as in Sec. 1.1.5. 



2.1 The glass transition of hard spheres in high dimension 

Hard spheres in the limit of large spatial dimensions provide an excellent opportunity for an analytic 
solution covering many aspects of the liquid and glass physics [139, 62, 140, 141]. The reason why this 
limit is solvable is geometric: consider three spheres A, B and C, with AB and BC in contact. What 
are the chances that A and C will themselves also be in contact? In high dimensions, vanishingly 
small. This led to the realization [139] that all terms in the virial expansion above the second could be 
neglected in high dimensions, as they involve geometrically heavily suppressed "coincidences", leaving 
one with only the two first terms of the series for the entropy 

S[p(x)} = J d d xp{x)[l - logp(x)] + l - J d d xd d y p(x)p(y)F(x - y) . (2.1) 

Here F(x) = e~ v ( x ) — 1 = — 9(a — \x\) is the Mayer function for e — > oo, i.e. for hard spheres, 
to which we restrict in this section. For the liquid phase, one has p(x) = p, a constant, and the 
above expression gives the liquid entropy and coincides with the Van der Waals equation. Non- 
uniform phases are described by solutions of the stationarity equation dS/dp(x) — that are not 
constant. Indeed, Eq. (2.1) provides the exact expression of the density functional for hard spheres 
when d — > oo, which has basically the form of Eq. (1.24). Actually, at very high densities Eq. (2.1) 
would have corrections, however these corrections are exponentially small in the density region that 
will be relevant in the following [141]. The liquid phase stays metastable at higher densities, where 
one expects the thermodynamics to be dominated by a modulated, crystalline phase, which is however 
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only known in small dimensions [142]. Despite these simplifications, the problem remains much more 
difficult than the study of simple mean-field models such as the p-spin model. The key difficulty is 
that, while for spin glass models the order parameter for the glass transition is the local Edwards- 
Anderson overlap, for particle systems the caging order parameter is a non-trivial function of space, 
the non-ergodic parameter defined in Eq. (1.5), or its Fourier transform, the long time limit of the 
van Hove function G s (r,t) [56]. 

2.1.1 Early theoretical results 

An early attempt to solve the amorphous DFT equations based on Eq. (2.1) was performed in [62]. 
First, a Gaussian density profile as in Eq. (1.25) was assumed. Next, it was assumed that the equilib- 
rium positions have the structure of a typical liquid configuration, hence So(q) — N" 1 J2j k e lq ^ Rj ~ Rk ^ ~ 
S(q) = 1/(1 — pF(q)) where F{q) is the Fourier transform of the Mayer function. Finally, a small p 
development was made to obtain So(q) <~ 1 + pF(q). Under these assumptions, a closed equation for 
the cage radius A was obtained. This equation has a finite solution when p > p d ~ 4.1 d2~ d 1 which 
defines the dynamical transition. In the same paper [62], the behavior of the MCT equations (1.27) in 
large d was considered. Under a Gaussian assumption for the non-ergodic parameter, f(q) ~ e~ Aq , it 
was shown that the MCT equations lead to the same equations for A of DFT. Hence, it was concluded 
that static and dynamical treatments lead to consistent results for d — > oo like in the p-spin model, 
and consistently with the RFOT theory. 

After the replica theory of [132] was developed, it was natural to apply it in larger dimension [143]. 
It was realized that a replicated version of Eq. (2.1), in which the coordinate x describing a single 
atom is replaced by a "molecular" coordinate x — {x\ ■ ■ -x m }, provides an exact description of the 
molecular liquid, for the same reasons why Eq. (2.1) provides an exact description for m = 1. The 
geometrical argument carries over to the molecules because when d — > oo the molecules become 
extremely compact, with a cage radius A oc 1/d 2 . Hence p{x) is replaced by p(x), with the advantage 
that the system of molecules is translationally invariant. A Gaussian assumption for p(x), of the form 



allows one to obtain an exact expression of the replicated entropy as a function of A. This procedure 
leads however to an equation for A that is different from the ones derived in [62] from DFT and MCT. 
In particular, it is found that tpd <~ 4.8d2~ d . A comparison of these equations seems to indicate that 
the RT equation contains more "diagrams" than the DFT and MCT ones. Indeed, an expansion of 
the RT equation in powers of F(q) gives back the DFT-MCT equation of [62] at the lowest order. 

A few years later, two groups independently found an asymptotically exact solution of the MCT 
equations (1.27) when d — > oo [104, 117]. The main features of this solution are that the f(q) is 
markedly non-Gaussian and it remains so when d — > oo, and the dynamical transition is found to 
scale as (fd oc d 2 2~ d , which is much larger than the previous results obtained from the Gaussian 
assumption. An inconsistency was pointed out in [117], namely that MCT gives a negative Van Hove 
function G(r, t) when d is large enough. 

At this point, one is faced with an unexpected puzzle: why in the limit d — > oo, where mean field 
theory should be exact, the three theories that are thought to be different facets of the unified RFOT 
theory give such different predictions? Looking carefully at the three theories, it is important to keep 
in mind that: 

• The DFT result is based -besides the Gaussian assumption- on a strong assumption, that 
So(q) ~ S(q). This is far from being obvious. In fact, when the pressure of the glass goes to 
infinity and A — > 0, the particles are stuck in their equilibrium positions Rj. But these have then 
to be strongly correlated, because of the requirement of mechanical stability of the packing. This 
will be further discussed in Sec. 2.2. In any case it is possible that the assumption So(q) <~ S(q) 
is too strong. 

• MCT has not yet been derived from a systematic large d expansion, therefore it is hard to tell 
whether the corrections to Eq. (1.27) become small when d — » oo. One could add terms to the 
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MCT equations that preserve the "universal" MCT phenomenology (i.e. the power-law scalings 
and the associated critical exponents) , but would change the "non-universal" predictions for the 
form of f(q) and the value of ipd [108]. At present, there is no reason to believe that the MCT 
equations (1.27) should be exact when d — > oo. 

• The RT results are obtained under one single assumption, that of the Gaussian form of the 
molecule distribution. They could therefore only fail if the Gaussian assumption were not 
verified. 

2.1.2 Numerical simulations 

In absence of any other theoretical indication of a way to resolve this inconsistency, one needs some 
numerical results to identify the correct path to follow. This was the main motivation to start a 
series of simulations trying to reach the largest possible dimensions in order to obtain indications on 
the scaling of the transition density and the shape of the non-ergodic parameter [144, 129]. Previous 
investigations had shown that crystallization is strongly suppressed when dimension is increased above 
3 and is basically impossible to observe on simulation time scales for any d > 4 [145, 146]. Hence, for 
these studies one can safely use monodisperse systems without having to worry about crystallization. 
Moreover, there was indication that the MCT predictions are very accurate in d = 4 [147]. 

Compression at constant rate [144] 

In a first set of simulations [144], we performed event-driven molecular dynamics simulations [148] 
of N identical hard spheres of diameter a, enclosed in a periodic box of volume V. A modified 
Lubachevsky-Stillinger algorithm was employed, which starts from a low-density gas (small ip) and 
density it by growing uniformly the particles at a constant rate a, reported here in standard reduced 
units such that the particle mass is m = 1, the temperature is ksT = 1, and a = 1. Time evolution 
stops when the system reaches a high reduced pressure p = /3P/p = 10 3 . Systems with N — 8000 
are simulated for d < 9 and larger ones for d = 10-13. These sizes ensure that even when the 
system is in its densest state the box edge remains larger than 2a, which prevents a particle from 
ever having two direct contacts with another one. There are strong reasons to believe that although 
relatively small these sizes nonetheless provide a reliable approximation of bulk behavior. First, with 
increasing d the box edge becomes less representative of the overall box size. The largest diagonals 
are \fd larger and there are many more diagonals than edges. Second, by analogy to spin systems, 
mean- field arguments indicate that for d — > oo, a hypercube of side two is sufficient to capture the 
full thermodynamic behavior. Even at the critical point, finite-size corrections are proportional to 
1/N S , where the exponent S is model dependent (e.g. 1/2 at the ferromagnetic transition), and do 
not directly involve the edge length L [50] . Similar results hold for dimensions greater than the upper 
critical dimension, where the exponents coincide with the mean field ones. Third, the fluid structure 
is expected to become uniform at ever smaller distances with increasing d [140, 141, 145]. Nearest- 
neighbor ordering should thus mainly be influenced by particles in contact or nearly so, with the rest of 
the fluid acting as a continuum. Indeed, in the fluid phase, finite volume corrections are proportional 
to the pair correlation h(L), and at fixed L, h(L) goes to zero exponentially with d [141]. The validity 
of these rationalizations has been satisfactorily tested by simulations in d = 8, where no finite size 
effect has been detected [144]. 

The compression results for d = 9 shown in Fig. 2.1 are representative of the behavior observed 
for all d > 3. The system first follows the equilibrium fluid equation of state (EOS) at low density 
and falls out of equilibrium at high density. Beyond this point, the pressure increases faster than in 
the equilibrium fluid and ultimately diverges at packing fraction <Pj(&). A Carnahan-Starling form 
PRuid(f) = 1 + 2 d ~ 1 </? captures well the pressure growth with if in the fluid regime (Fig. 2.1), 

provided, for each d, that one fits Ad to the data from the slowest compression rate available. Note 
that the contribution of the Ad term is found to decrease with increasing d, consistently with the 
asymptotic result of Eq. (2.1). In the high-density non-equilibrium regime, compaction runs with 
different & follow separate branches along which the pressure evolution is dominated by the expulsion 
of free volume [38]. Upon approaching jamming, the pressure is well approximated by the free volume 
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Figure 2.1: (Left, from [144]) Different compactions of N = 8000 particles in d — 9. With growing ip, 
the pressure first evolves like the fluid EOS then like a free volume EOS. In the inset, the fitted glass 
EOS parameters are reported as functions of a and the power-law extrapolation to a —¥ is shown. 
The resulting EOS for the threshold glass is reported in the main panel. It crosses the liquid EOS at 
(Pd, and its pressure diverges at <p t h- The prediction of replica theory for the ideal glass EOS is also 
reported for comparison: it leaves the liquid EOS at the Kauzmann transition tpx, and its pressure 
diverges at the glass close packing point ipccp, see Sec. 2.2. 

(Right, from [129]) Power-law fits (lines) to the vanishing diffusivity improve with <i=4-9 (different 
symbols), spanning over three decades of D. In the top inset, the resulting numerical tpd values, which 
coincide with the compression results from [144], are reported as triangles. The MCT prediction for 
ifd is reported with circles for comparison. In the bottom inset, the value of the exponent 7 from 
numerics (triangles) and from MCT (circles) is reported. 



form in which X/p is linear in tp: Pf v (a, ip) = rfyj ^ g ^_^ CT ^ . Both f(&) and tpj(&) are extracted from 
fitting the simulation data for p > p min ~ 50. 

These numerical results qualitatively agree with the RFOT scenario [1]. According to the theory, 
the compression curves for moderately small a should converge as a power law & a to a "threshold" 
glass. The dynamical transition density <pd separates the equilibrium fluid from this glass. Equilibrat- 
ing the liquid at ip > (pd requires much slower compaction rates, a < exp(— d), such that activated 
events allow the system to remain in equilibrium. In this regime one would expect a crossover to 
a logarithmic dependence of the glass EOS on <r, e.g., of the form 1 / 1 log(er) | , but it turns out that 
this regime is out of reach of present time numerical simulations. Indeed, the inset of left panel of 
Fig. 2.1 shows that the reported f(&) and p>j(&) are both linear as functions of & a with a ps 0.5. The 
exponent is expected to weakly depend on dimensionality, especially at low d, but this value is within 
the numerically reasonable range for all systems studied. Because the subsequent analysis is rather 
insensitive to the precise value of a, for simplicity it is kept constant. Extrapolating the results to 
(7 = gives the parameters iptt an d fth for the threshold glass free volume EOS, also reported in 
Fig. 2.1. Its intersection with the fluid EOS then provides a numerical estimate for tp^- These values 
are very close to previous numerical estimates in d — 3 [149] and d — 4 [147], which demonstrates 
the coherence of this analysis in low dimensions. The results for pd reported in Fig. 2.1 show that its 
dimensional evolution is smooth. Interestingly, dimensions where the crystal structure is singularly 
dense, such as d — 8 and d — 12 [142], do not present any echo of that singularity, which illustrates 
the smooth d dependence of the fluid structure and its independence from the crystal. 

Equilibrium dynamics in the liquid phase [129] 

Having obtained an estimate of tp d as a function of dimension, one can proceed to investigate the 
equilibrium dynamics in the liquid phase on approaching ipd- This has been done with the same 
algorithm as before, but stopping the compression at the desired ip < ipd, letting the system equilibrate, 
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Figure 2.2: (From [129]) The self Van Hove function, computed at the mid-plateau point tmp and at 
a density such that 2dD = 0.0008, hence close to ipd, for several dimensions d — 4,5,6,7,8. Dashed 
lines are the MCT predictions for d = 4, 6, 8. A fat non-Gaussian tail is clearly visible and is not 
reduced upon increasing dimension. 



and then performing equilibrium measurements of dynamical quantities. In particular, the average 
mean-square displacement (MSD) defined in Eq. (1.3) was obtained. At times shorter than the 
collision time, MSD displays a ballistic regime (Ar 2 (t)) = dt 2 , and at long times it has a diffusive 
regime (Ar 2 (t)) — 2dDt. Fitting these numerically determined D to the MCT power-law form 
D ~ (p — Pd) 1 is reasonably good for D < 0.005, and improves with increasing d. The resulting 
values of 7 and ipd are reported in Fig. 2.1. For d > 5, the full accessible dynamical range studied 
is used, spanning up to three D decades. The values of ifd obtained in this way up to d — 9 are 
consistent with the ones obtained from the constant compression data. The values of 7 are reported 
in Fig. 2.1. 

The cage form is described by the self part of the van Hove function 



which in the ballistic and diffusive regimes is well-approximated by a pure Gaussian. The logarithmic 
mid-point of the caging plateau tmp is chosen at the mid-time on a logarithmic scale, intermediate 
between the ballistic and the diffusive extrapolations of the MSD, and provides an estimate of the 
plateau regime discussed in Fig. 1.2 (see [129] for details). The results close to cpa, reported in Fig. 2.2, 
clearly shows that the cage is not Gaussian, with a fat tail at large r that does not shrink upon 
increasing dimension. By contrast, a comparison with the D5 crystal in d — 5 and the E8 crystal in 
d = 8 (unpublished) shows that the cage is there very well described by a Gaussian form. 

Discussion 

In summary, numerical results indicate that: 

• The MCT qualitative predictions are very well obeyed by the data, as previously observed in 
d = 4 [147]. The diffusion constant vanishes as (p — pd) 1 ■ The vanishing of the diffusion 
is accompanied by a growing caging plateau. It remains to be checked whether the approach 
and departure from the plateau define power-law regimes in time with exponents a and b, and 
whether the relation between 0,6,7 is consistent with MCT predictions. 
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• Dynamical heterogeneity is reduced with increasing dimension [129]. Together with the previous 
observation, this suggests that indeed when d increases the system approaches a mean field 
regime well described by the RFOT scenario. 

• Yet, the quantitative MCT predictions for 7 and the cage radius become increasingly poor 
with increasing d. This is not a problem for the RFOT scenario: it is reasonable to expect that 
a different theory with a similar structure to MCT but a different kernel exist. Such a theory 
would provide the same qualitative predictions as MCT [108], and exact quantitative results for 
ifd, 7 and the cage structure. But this theory has yet to be found. 

• The Gaussian assumption for the cage form done by RT and DFT is invalidated by the data, 
calling for a revision of these theories. Unfortunately, the predictions for tp^ of RT in low d 
are not reliable because of the inaccuracy of the small cage expansion (the replicated HNC 
calculation has not been performed yet). However, the result of RT for px should provide an 
upper bound for p d , and this is consistent with the data. Also the predictions of RT for the 
ideal glass EOS (Fig. 2.1) are consistent with the data, as the ideal glass EOS should be similar 
to the threshold glass EOS but shifted in density. 

2.1.3 New theoretical results 

The numerical results obtained in [144, 129] give enough suggestions on how to proceed to improve 
the theory. 

On the side of replica theory, a first step has been obtained in [150], where the exact description 
of the molecular liquid in the high-dimensional limit, without further assumptions, was found. The 
strategy is the following: a generic replica problem is written in terms of the order parameter p(x). 
The x a are vectors in d-dimensional space, and one is looking for a solution that is statistically 
translationally invariant and isotropic. The only possibility with these properties is that p(x) depends 
exclusively on the \x a — X(,\ 2 = q aa + q^ — 2q a { )1 where q a b = x a -Xb. All d-dimensional integrals may be 
expressed as low dimensional integrals in terms of the q a b, with volume factors, a simple generalization 
of spherical coordinates. Then, the dimension d appears explicitly, and the limit of large d may be 
taken in a straightforward way, by saddle point evaluation of the integrals. The result is that the 
Gaussian ansatz turns out to be, in a sense, exact: it gives for large dimensions the exact values 
for the thermodynamic quantities. The reason is that in the large d limit all the saddle points are 
dominated by a typical distance between atoms in a molecule, and therefore the exact shape of the 
cage function is irrelevant. Any reasonable ansatz would give the same results. This explains why 
the Gaussian results for the glass EOS, as well as the transition densities, are consistent with the 
numerical data, as shown in Fig. 2.1. The next step in this direction is to work out the prediction of 
RT for the cage shape and compare it with the results in Fig. 2.2. This work is currently in progress. 

On the side of dynamics, it might be possible to obtain an exact solution starting from the approach 
of [151] and following the same lines of the replica calculation. Also this work is currently in progress. 

In summary, although the physics in large d seems quite well understood, its mathematical de- 
scription is surprisingly complex. Yet a lot of progress has been made and it seems that we are on the 
right path towards a complete RFOT theory in this limit. 

2.2 Mean field theory of the jamming transition 

The main features of the jamming transition have been introduced in Sec. 1.1.5. Consider the simple 
model of harmonic spheres described by Eq. (1.11). The two control parameters are the packing 
fraction ip and the ratio of thermal energy to potential energy T/e (we always use k B = 1). As 
discussed in Sec. 1.1.5, the jamming transition happens out of equilibrium at a packing fraction ipj, 
when the system is inside its glassy phase. At this transition, the origin of the rigidity of the glass 
changes in nature [40]: below ipj, particles do not overlap and rigidity is entropic in nature, in such a 
way that pressure and elastic moduli are proportional to T; above tpj, particles overlap and rigidity is 
mechanical in nature, hence pressure and elastic moduli have a finite limit at T — > 0. Hard spheres have 
e — > 00, hence T/e — > 0, and they lie on the horizontal axis in the phase diagram of Fig. 1.6, but they 



2.2. MEAN FIELD THEORY OF THE JAMMING TRANSITION 



31 



exist only for ip < ipj. At ipj, the reduced pressure of hard spheres, defined by p = P/(pT), diverges, 
see Fig. 2.1. This can be understood in two equivalent ways. From the hard sphere perspective, ipj 
is the point where the spheres get in contact with their neighbors. The cage radius goes to zero, and 
collisions with the neighbors become infinitely frequent. If T is fixed, the momentum transfer per 
unit time between neighbors, which is proportional to p, diverges. From the soft sphere perspective 
(i.e. at finite e), the pressure becomes non-zero at ipj even at T = 0, hence p oc P/T = oo in the 
whole region above ipj. In other words, ipj is identified with the close packing point [152, 153], the 
maximal density at which the spheres can be packed avoiding overlaps. This has been called random 
close packing, <prcp, in the literature. In present notations, ipj and fpcp coincide. 

Independently from these motivations originating in the physics of amorphous materials, the prob- 
lem of packing spheres in d dimensions is intimately related to several important mathematical prob- 
lems, notably in the context of signal digitalization and of error correcting codes, especially when d 
is large [142]. Although the densest packings are known (either rigorously or at the level of good 
conjectures) in low dimensions, as soon as d becomes large, almost no results are available, see [142] 
and the very nice review paper [154]. Providing a "physical" estimate of the closest packing den- 
sity in large dimensions would be very interesting for the mathematical community. Attempts in 
this direction have been performed, and even if intuition suggests that lattice structures should fill 
space better than amorphous ones, some have suggested that for large d disordered packings may 
be at least as efficient as lattice-based versions [155]. Although the density of lattice geometries is 
strongly d dependent [142], the fluid structure and its properties are much less sensitive to d, once 
d > 3 [146, 144, 129]. A general understanding of the jamming transition in arbitrary d would thus 
also clarify the density scaling of amorphous packings and their potential mathematical role. 

In summary, the jamming transition is interesting for wildly different problems, ranging from the 
physics of foams, pastes, emulsions, powders and other granulars, to the physics of thermal glasses 
in the limit of low temperature, and to mathematics and information theory [1, 156]. Despite all 
this interest, and despite the fact that the random close packing concept has been proposed long 
ago [152, 153], a first principle theory of jamming has not been developed. "First principle" here 
means a theory that is able to predict the existence of the transition, its location in density, and 
the associated structural changes, from the sole knowledge of the particles' interaction potential. A 
first step in this direction was obtained -unexpectedly- as a side product of the replica theory of 
hard sphere glasses in [132]. Subsequent improvements have been made in [1, 2]. This theory will be 
reviewed in the rest of this section, highlighting its successes and the many problems that are still 
unresolved. 

2.2.1 Is random close packing well defined? 
Predictions of replica theory 

How does the jamming transition appear in the replica theory? Recall from the discussion of Sec. 1.2.2 
that one has to introduce to replicas as a tool to take into account the existence of multiple amorphous 
glassy states, m is a free parameter and in the end one has to optimize the free energy with respect 
to it, the optimal value for a given state point (<p,T) being called m*(<p,T). It is expected that 
to* = 1 in the liquid phase, and m* continuously becomes smaller than 1 at the glass transition [87]. 
What came as a surprise in [132] is that, for the hard sphere system, m*(ip,T = 0) continuously 
decreases starting from 1 at (f = (fx, and eventually vanishes linearly at some "(ideal) glass close 
packing" (GCP) density (ficcp, m * ( ( Pgcp,T = 0) ~ \<p — (pccp\- This is crucial, because within 
replica theory the reduced pressure p of the (ideal) glass is proportional to 1 /to* , hence the vanishing 
of to* drives a divergence of p and signals the jamming transition! Note that this result was not put 
by hand in the theory: the calculation of the replicated free energy can be done starting from the 
interaction potential, as a function of (ip, T) and to, and it is the process of optimization with respect 
to to and the vanishing of to* that is behind the appearance of the transition. A posteriori, the 
fact that the jamming transition is related to the vanishing of to is obvious thanks to the connection 
with combinatorial optimization problems where a similar transition happens (and is called the SAT- 
UNSAT transition in that context [157]). This connection was explored in details in [158, 159]. Note 
also that the lattice glass model of [65, 160] is a perfect example of this connection. 
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However, what is the relation between <pgcp ancl the random close packing of Bernal [152, 153]? 
The GCP is the point where the pressure of the ideal glass diverges. However, as discussed in Sec. 1.1.1, 
even if one accepts the mean field prediction that such ideal glass really exists, producing it in a 
reasonable amount of time is impossible both in numerical simulations and experiments. The liquid 
will always fall out of equilibrium much before ipx is reached, and a non-equilibrium glass of lower 
density will be obtained. This is well illustrated in Fig. 2.1. Hence, in real- life experiments and 
numerical simulation, one will be stuck in out of equilibrium states, and in general the random 
close packing point obtained in a given experiment will be at lower density, ipj = ifiRCP < fGCP, 
corresponding to a system falling out of equilibrium at some (p g < ipx- 

It might seem that the replica description of the GCP is therefore useless. Fortunately, this is not 
the case. Indeed, as discussed in Sec. 1.2.2, RFOT theory predicts the existence of (exponentially) 
many metastable states, each corresponding to a different glass, with different free energies. And 
replica theory allows one to compute their properties. Performing the computation [1], one discovers 
that each different glass will jam at a different tpj. Therefore, replica theory predicts the existence 
of a range of densities ipj e [</?™ m , feep] m which amorphous jammed packings exist. Each of them 
corresponds to jamming of a different non-equilibrium glass. Therefore, replica theory predicts that 
the jamming (or RCP) point is not well defined: it depends on the protocol that one uses in the 
experiment. Actually, the whole phase diagram in Fig. 1.6 (i.e. the transition line T g (<p) and the 
jamming point ipj) depends on the particular metastable state (as parametrized for instance by the 
equilibrium positions Rj) that is selected in the experiment under consideration. So what is left of 
the predictive power of replica theory, besides setting upper and lower bounds for ipj? Here comes 
another crucial observation: replica theory predicts that all amorphous packings that correspond to 
jamming of a metastable glass have equivalent structural properties [1]. This means that even if the 
jamming density ipj is protocol-dependent, the structural signatures of jamming are universal. 

This prediction of replica theory should not be confused with the idea that RCP is not well defined 
because one can construct packings with arbitrary degree of crystallization [161, 162]. Although this is 
certainly true, the presence of such microcrystallites is detectable using structural measures (so-called 
order metrics [161]). The prediction of RT is that there is an interval of perfectly amorphous packings. 
Of course those can be mixed with crystalline regions to produce any intermediate density up to the 
crystal close packing. 

Numerical tests 

This prediction has been successfully tested in a number of independent studies [35, 164, 165]. How- 
ever, these studies were limited to d = 3, hence subject to the criticism that crystallization might play 
a hidden role in the density variation [161]. Moreover, the density range of amorphous packings in 
d = 3 is quite small [164]. Hence, more recently [163], a systematic study as a function of dimension 
was performed, as part of the previously mentioned effort towards a systematic understanding of the 
RFOT theory in large d [144, 129]. Amorphous packings were prepared using two different numerical 
protocols. 

(i) Approaching jamming from densities below it by Lubachevsky-Stillinger (LS) compressions of 
hard spheres undergoing Newtonian dynamics while a grows at a fixed rate &, as previously 
discussed. The compression stops when particles are in contact, defining the packing fraction 
ifp, at which the reduced hard spheres pressure becomes infinite. 

(ii) Approaching jamming from densities above or below it by minimizing the energy E(X) = 
J2i<j v (\ x i ~ x j\) °f a S0I t sphere system with interaction potential v(r) in Eq. (1.11) [31]. One 
starts from random sphere configurations X + = X_ and uses ct_ and a + that bracket the 
jamming point, in such a way that E{X_) = at cr_, while E(X + ) > at a + . The onset of 
non-zero energy is found using an iterative bisection method. At each iteration, an intermediate 
value <7 m is chosen, and the energy at a m is minimized starting from either X + or X_. This 
choice defines the procedure to approach jamming from above or from below, respectively. Calling 
X m the configuration obtained after minimization, one substitutes a + = er m and X + = X m if 
E(X m ) > 0, while er_ = er m and A"_ = X m if E(X m ) = 0. The procedure stops when E(X + ) 
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Figure 2.3: (From [163]) Main panel: the extrapolated jamming density <fp~^° following the protocol 
described in Fig. 2.1 (full line and crosses), and compared with the prediction of RT for (ficcp (dashed 
line). Top inset: the range of jamming densities (squares) is compared to (^™ ax (circles) and v?™ m 
(triangles). Note that tp™ ax ~ ipp =3xlQ and </?™ m ~ (fp =3xl ° . Bottom inset: the increase of ip\ 
with cr + , in terms of the initial effective packing fraction, in d — 3. 



is below a given threshold (of the order of numerical accuracy), and we call the final density 
\ in the procedure from above, and (p\ in the procedure from below. From above, the energy 
vanishes with e = E/N ~ Aip 2 and the static pressure P ~ Aip, where Aip is the distance from 
jamming [31] 

For protocol (ii), both <j± have no measurable effect on ip\. One can formally define <pf m = 
min^ ip\{a±), but any reasonable <7± results in the same final density. By contrast, <p\ strongly de- 
pends on er + (Fig. 2.3), but is also independent of <r_. One can therefore define <^™ ax = max CT± ip^(a±). 
A practical way of constructing both ip™ ln and ^™ ax is to run the energy minimization (respectively 
from below and from above) starting from er_ = and a <r + large enough that ip* saturates to its 
maximum; intermediate packing fractions can then be obtained by reducing a + (Fig. 2.3). It is found 
that [<p™ n , </?™ ax ] roughly corresponds to {ipp~ , ipp + ] with cr_ ss 3 x 1CT 2 and & + « 3 x 1CT 4 , and 
that large a generate mechanically unstable packings. In other words, although the jamming density 
sensitively depends on the protocol's specific parameters, the jamming density range is much less 
sensitive to the protocol's class. The resulting density range is found to grow steadily from about 
2% in d = 3 to nearly 10% in d = 11. Therefore, the similar observations made for d = 3 binary 
mixtures [35, 164, 165] are confirmed, but in larger dimension the effect is much larger and leaves no 
room for ambiguity, also because crystallization is totally absent as soon as d > 3 [146] . 

The value of <^™ ln ~ ^p^° are both lower than (fccp predicted by replica theory (Fig. 2.3). This 
is expected because, according to the analysis of the glass problem (Sec. 1.1.1) and of analogous 
combinatorial optimization problems [159, 158], achieving a packing density of fccp requires at least 
to equilibrate the liquid close to (fx, which can only happen through activated processes and can be 
expected to require a time increasing exponentially with d (or with N in mean field models), hence 
being out of reach of numerical simulations. Unfortunately, a reliable prediction from RT for the 
minimal value of the jamming density cannot be obtained, for similar reasons as for cp^. 

2.2.2 Finite temperature scaling relations 

As shown in Fig. 1.6, the jamming transition exists only at strictly zero temperature (here we consider 
soft spheres with a finite e). In fact, as soon as temperature is non-zero, overlaps between particles 
are allowed by thermal fluctuations, and a strict distinction between zero energy non-overlapping 
configurations and positive energy overlapping configurations is not possible. However, if temperature 
is low enough, the proximity of the jamming transition induces scaling relations for all thermodynamic 
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Figure 2.4: (From [2]) Results for soft spheres in d = 3. (Left) Optimal number of replicas to* as 
a function of T, for several volume fractions. The inset shows the T — > limit of to* as a function 
of the volume fraction. (Right) Scaling of m* (T, tp) /VT as a function of the rescaled temperature 
I <P — <Pgcp I /VT for tp = 0.58, 0.59, . . . , 0.70. The asymptotic forms corresponding to Eq. (2.4) are also 
reported. 



quantities. Here these scaling relations are discussed. For this discussion, one can restrict to a given 
non-equilibrium glass with given jamming density tpj (for instance, the ideal glass with jamming 
density tpccp) because the results are independent of the choice of glass state and the corresponding 
tpj. 

It turns out from the solution of the replica equations [2] that all the thermodynamic quantities 
are singular at tpj and T — 0; for instance m* is finite below tpj while it vanishes proportionally to T 
above tpj . The reduced pressure is also finite below tpj , while it diverges at tpj and is formally infinite 
above tpj. Both the energy and the pressure vanish below tpj while they are finite above tpj. From 
this observation, and because all quantities are analytic at finite T, it follows that all these quantities 
must satisfy scaling relations if T is small enough and tp is close enough to tpj [166, 34, 35]. In the 
replica theory, all the scalings are driven by the one of m*(tp,T), for which one finds 

rh_(x — > oo ) = Jlx , 

(2.4) 

fh + (x — > oo) = — . 

TX 

This ensures that in the hard sphere limit T — > and tp < tpj, one recovers m* ~ \tp — tpj\, while for 
tp > tpj, one finds that to* = T/\tp — tpj\. Note that p, and f can be computed explicitly from the 
theory and they slightly depend on tpj [2]. This scaling is shown in Fig. 2.4. 

Qualitatively similar scaling forms (in some cases with different exponents) apply to thermody- 
namic quantities such as the energy, the pressure, or the complexity. Although the corresponding 
scaling plots are not shown explicitly, the energy and the pressure of the glass are plotted and com- 
pared with numerical data in Fig. 2.5. The latter have been obtained by preparing an amorphous 
packing through the energy minimization protocol discussed above, and then performing standard 
Monte Carlo simulations at finite temperature in the vicinity of tpj ~ 0.643, and checking care- 
fully that no structural changes are observed, because in the simulation one is always in the regime 
T <C Tg(tp) [2]. The agreement between theory and numerical data, with a rescaling of dtp. discussed 
in [2] that is needed to correct for a thermodynamic inconsistency of the theory, is nearly perfect. The 
scaling around jamming is clearly visible in the figures. For instance, the average energy U g i ass (T, tp) 
tends to a finite value for Stp > 0, while for Stp < it goes to zero as a power law, because in this case 
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Figure 2.6: (From [2]) Left panel: evolution of the maximum of the glass pair correlation function 
with T and ip. Right panel: average number of contacts z as a function of distance from jamming Sip, 
compared with numerical data (as in Fig. 2.5), for several temperatures. 



the system becomes a hard sphere glass. Similarly, the reduced pressure is finite for Sip < 0, while it 
diverges proportionally to f3 for Sip > 0, because in this case the pressure is finite at zero temperature. 
At finite temperature, the curves interpolate between the two regimes. Similar scalings have been 
discussed in [167]. 

2.2.3 Scaling relation for the structure 

At the jamming transition, a marked structural change happens, as illustrated in Fig. 1.6. Particles 
are not in contact below ipj, they just touch at ipj, and they overlap above ipj. This is clearly revealed 
in the pair correlation function g(r) = (pN) ' \%2i^jS(r + Xi — Xj)) -the Fourier transform of the 
structure factor S{k)-, and in the cumulative structure function 

Z(r) =pil d f dss^gis) , (2.5) 
Jo 

where is the e?-dimensional solid angle. The function g(r) is proportional to the average number of 
particles at distance r from a given particle, while Z(r) counts the average number of particle centers 
in a shell of radius r around a given particle. 
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Figure 2.7: (From [2]) Pair correlation g(r) in the glass phase just below (left panel) and just above 
(right panel) jamming, predicted by theory (full lines) and measured in numerical simulations (sym- 
bols). 



Number of contacts 

In fact, g(r) develops a delta peak in r = a when T = and <p = tpj, corresponding to the touching 
neighbors. The coefficient of the delta peak, which is also equal to ~z = Z(r = a), gives the average 
number of neighbors. The maximum of g(r) as well as z are reported in Fig. 2.6 as a function of 
(T, ip) [2] . It is indeed observed in the theory that g m ax diverges on both sides of the transition at 
T = as g max ~ \<p — <fj\ , while this divergence becomes a smooth maximum at finite T near the 
transition, whose position shifts with temperature. This behavior compares very well with numerical 
[38, 39] and experimental [36, 168] observations. 

Similarly, z is a smooth function of density at T > 0, but it develops a jump at T — and 
ip = Lfj. The limit z(T = 0, ip — > <p~j) = 2d corresponds to the so-called isostatic condition. This is a 
crucial property of jammed amorphous packings [31, 169, 33, 32]. Indeed, it is known since the work 
of Maxwell that isostatic systems are at the verge of mechanical instability: removing one contact 
creates one "soft" mode that allows particles to slide without paying any energy cost. This is believed 
to be the origin of all the peculiar mechanical properties of nearly jammed systems [31, 170, 33, 32]. 
It is therefore quite surprising that a mean field theory, that is a priori completely unrelated to the 
detailed mechanical structure of the system, is able to predict that z — 2d at the transition. However, 
an important discrepancy between theory and numerical data is revealed in Fig. 2.6. In fact, while 
the theory predicts z — 2d oc Sip for T — and 5<p > 0, numerical data show that z — 2d oc y/5(p (we 
will discuss this point later). 

Scaling of the pair correlation around the delta peak 

The shape of g(r) is deeply modified around jamming with respect to the liquid phase, due to the 
development of the contact structure. The most spectacular feature is the emergence of the delta peak 
at T — ¥ and tp — >• ipj. Replica theory gives a prediction for the shape of the delta peak in the whole 
region for T > and (p ~ cpj, which is tested against numerical data in Fig. 2.7 and is shown to work 
quite well [2]. 

A more detailed investigation can be performed at T = 0, approaching ipj from above and from 
below using protocols (i) and (ii) described in Sec. 2.2.1 [163]. For both protocols, at fixed T = 
and (p ~ ipj with Aip = \5<p\ — \<p — <pj\, it is found that Z(r) jumps from to the isostatic average 
number of contacts z ~ 2d around r = a, on a scale r — a cx Aip. This regime corresponds to the 
contact delta peak of g(r). The plateau is extended by a second power-law regime that corresponds to 
particles in "quasi-contact", carrying no force at jamming. In this regime the scaling is the same for 
both protocols and Z(r) displays the same growth until it reaches the trivial large r regime. A third 
intermediate regime is required to match the first two around the plateau, as illustrated in Fig. 2.8. 

Consistency sets clear scaling requirements for the different regimes. When approaching jamming 
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Figure 2.8: (From [163]) Schematic of Z(r) when approaching jamming (b) from below and (a) from 
above. Three distinct scaling regimes can be identified. The first regime is related to the growth of 
Z(r ) from to z. It corresponds to interparticle gaps h = \r — a\ ~ Atp, and hence to particles that are 
in contact when Aip — > 0. The last regime corresponds to gaps h that remain finite for Aip — > 0, and 
therefore these particles remain separated at jamming. These small gaps, Z(r) — z oc h 1 ^ 01 , nonetheless 
result in a "quasi-contact" regime. The intermediate regime corresponds to gaps h ~ Aip^. It matches 
the two other regimes and disappears when Aip — > 0. 



with protocol (i), Z{r) = for r < cr, and for r > cr, the pressure p oc Atp~ l parametrizes the scaling 
function for Z(r > cr). A first scaling regime r — a ~ sees Z(r) grow from to the average number 
of "contacts" z as 

Z(r) = zZ_ [(r - cr)p/cr] (2.6) 

with Z_(x) ~ 1 — Cx^ 1 ^ 6 when x — > oo for a constant C [38]. Force-bearing contacts are only 
observed at jamming proper, but their signature develops asymptotically. A second regime for finite 
r — a has 

Z(r) = z + C'(r-a) 1 - a , (2.7) 

where C is a constant. At jamming, these nearly touching "quasi-contacts" carry no force. For large 
r, a trivial regime develops independently of \Aip\. Matching the first two scaling regimes implies the 
existence of an additional intermediate regime for r — a ~ where 

Z(r) = z + p v -^n-[(r-a)p^/a] (2.8) 

with n < 1 and v < fi (Fig. 2.8). Consistency then requires that T-L^(x — > 0) oc — x^ 1 ^ and 
%_{x — > oo) oc with scaling relations v = afi and fx = (1 + 0)/(2 + — a). 

When approaching jamming with protocol (ii) from above, the remaining overlaps provide a scaling 
variable e — \ fe~dj (2e) oc |A<^| [2], In spite of the very different preparation protocol, similar structural 
regimes are identified. In the r < a contact regime, Z[r) grows from to I by a universal scaling 
function 

Z(r) = zZ + [(cr-r) £ -7cr] (2.9) 
with Z + (x) ~ 1 — C"x 1+e when x — > for a constant C" . For r > cr, here again 

Z(r) = z + C'(r-a) 1 - a , (2.10) 

hence the two regimes must be matched by an intermediate scaling function % + for r — cr ~ 

Z{r) = z + e^-' / U + [{r-a)E-»/(T] (2.11) 
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Figure 2.9: (From [163]) Left panel: scaling of Z(r) upon approaching jamming by hard sphere 
compressions (p oc |A(^| _1 — > oo) in d = 3. (a) For different p, Z(r) grows on a scale r — a ~ p^ 1 
to a plateau at z ~ 6, the isostatic value, (b) The small r — a regime shows the "contacts" scaling 
function Z-(x), which agrees well with the G-RT prediction (red line) at small r — a. (c) Rescaling 
Z(r) using /i = (1 + 0)/(2 + 9 — a) and v — ap highlights the behavior of the scaling function 
\H-(x) - H-(l)\ ~ lOx (brown line) along with the 9 = 0.30(3) (red line) and a = 0.40(1) (blue line) 
power-law regimes, (d) The "quasi-contact" regime with a = 0.40(1) (blue line) eventually gives way 
to the trivial r d regime as r ^> a [38] . 

Right panel: scaling of Z(r) upon approaching jamming by soft sphere energy minimization from 
above in d = 3. Recall that the energy e ~ |A0| 2 and the scaling variable e — yJe d/(2e) oc |A<^|. (a) 
For diminishing e, the height of the plateau converges to the isostatic value Z(a)=2d. (b) The small 
r < a regime shows the "contact" scaling function Z + (x), which agrees well with the G-RT prediction 
(red line) at small r — a. (c) Rescaling Z(r) using /j, = (1 + ff)/{2 + 9 — a) and v — a,[i highlights 
the behavior of the scaling function \K-ix) — H-(l) \ ~ 1.2a; (brown line) along with the 9 — 0.42(2) 
(red line) and the a = 0.39(1) (blue line) power-law regimes, (d) The "quasi-contact" regime with 
a ~ 0.39(1) (blue line) eventually gives way to the trivial r d regime as r > <r. 



with fj, > 1 and v < 1 (Fig. 2.8). Consistency here requires that Ti+(x —> — oo) oc — and 
H+(x — > oo) oc x 1 ~ a with v = ap and /i = (1 + 9) /{a + 9). 

These scalings are verified using numerical data from [163] in Fig. 2.9. The Gaussian replica theory 
provides precise and simple predictions for the contact regime scaling function Z±{x) [1, 2, 163]: 



Z-(x) = l-e^ x 
Z+(x) = 1 — erf(x) 



1 — erf 



^7T 



(2.12) 



These predictions are found to be extremely accurate when x is of order 1, but they fail to capture 
the power-law tails at large or small x (Fig. 2.9). Gaussian RT indeed predicts an exponent 9 = 
for both protocols, and completely misses the power-law divergence related to a, predicting a = 0. 
Note that, although the exponent 9 plays an important role as it will be shown below, the deviation 
between the RT prediction and numerical data is barely visible in Fig. 2.9, and it is not visible at all 
in Fig. 2.7 because it is hidden in the tails of the delta peak, which are also smoothed by thermal 
fluctuations. 

The origin of these discrepancies is not clear yet. One possibility is that they are due to the 
Gaussian assumption for the cage form, which was found to be erroneous as detailed in Sec. 2.1.2 
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Figure 2.10: (From [163]) Cumulative force distribution G(f) compared with the Gaussian RT pre- 
diction (solid line) (a) for HS with power-law tail exponent 9 = 0.30(3) (dashed line) and (b) for SS 
with power-law tail exponent 8 = 0.42(2) (dashed line). Inset: test of Eq. (2.13) in d=3. The left 
hand side (from distances) is plotted with points, the right hand side (from forces) is plotted with a 
solid line. 



and in [129]. It is possible that the non-Gaussian replica theory of [150] will lead to an improved 
prediction, however this has still to be worked out. The other possibility is that these discrepancies 
have their origin in the physics of the low frequency vibrational modes (often called "soft modes" ) 
that appear at the jamming transition due to the approaching mechanical instability [31, 170, 33, 32]. 
It is worth stressing that the analysis of Fig. 2.9 has been repeated in [129] in dimensions d = 3- • • 10 
and no difference has been found. The fact that these scalings persist up to high dimensions strongly 
suggest that the mean field theory should be able to capture them. However, how to include the soft 
modes physics in the replica theory is not clear yet. 

Force distributions and mechanical stability 

The scaling functions Z± (x) characterize the contact regime and are related to corresponding scaling 
functions for the probability distribution of inter-particle forces /. Here again, it is convenient to 
consider the cumulative distribution G(f) rather than the distribution itself. 

For hard spheres approaching jamming from below, the average force / cx P and the cumulative 
force distribution G(f) approaches a scaling function defined by G{yf) — > G-{y)- In the contact 
regime the force and distance distributions are also related through a Laplace transform 



as suggested in [38] and tested in Fig. 2.10. It follows that if Z_{x) ~ 1 — Cx^ 1 ^ 6 for x — > oo, 
then G-(y) ~ y 1+e for y — > 0. The low-force distribution is indeed consistent with G(f) oc / 1+ and 
9 = 0.30(3) (Fig. 2.10). For soft spheres approaching jamming from above, the pair potential sets the 
relation between the force and the pair distributions: it gives / = 2(e/cr)(l — r/a) for < r < a and 
zero otherwise, so G(f) = 1 — Z ^ a ~zlJ) t " 2t ^ ■ ^ n tne j ammm g limit G(2yee/a) — > G+{y) = 1 — Z + (y), 
and therefore G+(y) ~ y 1+e , as in the previous case. Indeed, here again, the low- force tail is consistent 
with 9 = 0.42(2) (Fig. 2.10). For both protocols, however, the regime intermediate between contacts 
and quasi-contacts results in deviations from this power-law decay at very weak forces away from 
jamming. 

The exponents 9 and a defined above might have an important role in the mechanical stability of 
the amorphous packings. In fact, it has been proposed by Wyart [171] that a > 1/(2 + 9) is required 
for mechanical stability. Both the soft spheres values (a = 0.39(1), 9 = 0.42(2)) and the hard sphere 
ones (a = 0.40(1), 8 — 0.30(3)), however, indicate a slight violation of this condition. A generalized 
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Figure 2.11: (From [172]) Packing fraction ipj as a function of 77 = 1/(1 + xr 3 ) at fixed r. Full 
symbols are numerical data from [172]. Open symbols are experimental results from [179]. Lines are 
predictions from theory. Note that the large r-small r\ region cannot be explored, because for such 
very asymmetric mixtures the large spheres form a rigid structure while small spheres are able to 
move through the pores and are not jammed [177, 178]. 



stability condition of the form a > (1 — S/2)/(2+6 — 8/2) [171] is consistent with the numerical findings 
for 5 > 0.2, but a direct test of this extended relation goes beyond the scope of the current analysis 
because the measurement of the value of 5 requires a completely different numerical procedure. 

2.2.4 Binary mixtures 

The replica theory can be extended to binary mixtures [172]. This is important because in these 
systems crystallization can be easily avoided also in d = 3; this allows one to study their glass transition 
in great detail, see e.g. [173, 174, 175, 122, 135, 34, 35, 165]. Jamming of binary mixtures has also 
been extensively investigated, see e.g. [176, 177, 178, 31]. Their investigation allows one to test the 
prediction of the theory concerning the variation with composition of packing density, structure, etc. 
Here results are reported for a binary mixture of two types of three-dimensional spheres p, = A, B in 
a volume V , with different diameter <r M and density p M = N^/V. Define r = <ja/o~b > 1 the diameter 
ratio and x = Na/Nb the concentration ratio; tp — paVz(cfa) + PbV$((Jb) the packing fraction; 
T] = pbVz(<jb)/v — 1/(1 + xr3 ) the volume fraction of the small (B) component. The equation of 
state used in [172] to describe the liquid is a generalization of the Carnahan-Starling equation. The 
theoretical results were compared with numerical results obtained from the Lubachevsky-Stillingcr 
algorithm discussed in section 2.1.2. 

In Fig. 2.11, the jamming density tpj is reported for different mixtures, putting together numerical 
results [172], experimental data [179], and the theoretical results. The latter have been obtained 
by fixing the complexity to a single value for all mixtures, such that for the monodisperse system 
ifj ~ 0.64. Note that a single "fitting" parameter allows one to describe different sets of independent 
numerical and experimental data. The prediction of the theory are qualitatively similar to previous 
ones [177, 178], but the quantitative agreement is much better. Interestingly, a similar qualitative 
behavior for ifd has been predicted by Mode-Coupling theory [173, 174]; although there is no a priori 
reason why the variation with mixture composition of ipa and <pj should be related, it is reasonable 
to expect that they show similar trends [174]. 

The average coordination numbers at ifj are denoted z^ u ((pj), and it was checked in [172] that 
their variations with tpj are negligible. They are reported in figure 2.12 for different mixtures. The 
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Figure 2.12: (From [172]) Partial average coordination numbers (small-small, small-large, large-small, 
large-large) as a function of volume fraction of the small particles r\ = 1/(1 + xr 3 ) for different values 
of r. Full symbols are numerical data from [172]. Open symbols are experimental data from [180]. 
Note that in the lower right panel a different scale is used for z; s . 



numerical values have been obtained in [172]. Experimental data from [180] are also reported in the 
right panel of figure 2.12. The total coordination is close to the isostatic value z = 6, which is the 
value predicted by the theory also for binary mixtures [172]. As it can be seen from figure 2.12, the 
computed values agree very well with the outcome of the numerical simulation, at least for r not too 
large, while some discrepancies are observed in the contacts of the large particles for large r. 

2.2.5 Discussion 

In summary, replica theory provides a quite satisfactory mean field theory of the jamming transition. 
First of all, the very fact that the theory predicts the existence of the transition from first principles 
is non trivial. Furthermore, the basic phenomenology of the transition, with the associated diverging 
peak of g(r) and jump of the number of contacts to the isostatic value, are reproduced. And finally, 
quite accurate predictions for the scaling functions Z±(x) can be obtained, and compare quite well 
with numerical data. 

Despite these successes, the level of accuracy of numerical simulations around the jamming tran- 
sition is so high that very detailed measurements can be made, and these reveal several discrepancies 
with the theory. In fact, the exponents 9 and a are predicted by the theory to be = a = which is 
wrong. Furthermore, the growth of the number of contacts above the isostatic limit is predicted to be 
linear, while numerics indicate a power-law growth with exponent close to 0.5. These discrepancies 
are crucial, because these non-trivial exponents have been related to the mechanical stability of the 
packings [170, 171], and to the proliferation of soft modes that is connected to the isostatic nature of 
the force network. 

The results of [163] indicate that all the non-trivial exponents are independent of dimension, 
therefore the nature of the jamming transition in d = 13 is exactly the same as in d = 3. This strongly 
suggests a mean field description and indeed it has been suggested [181] that the "upper critical 
dimension" for the jamming transition is 2 (although an equilibrium glass transition probably does 
not exist in d = 2 [135, 122]). Moreover, it has been shown [66] that a model of hard spheres on the 
Bethe lattice, which is therefore rigorously mean field in nature, reproduces all the phenomenology of 
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the jamming transition, including the soft modes physics. It seems therefore that one can construct a 
refined mean field theory, capable of predicting the non-trivial exponents and the associated peculiar 
mechanical features of nearly jammed packings. 

How can one do this in practice? Some different paths can be devised. The first one is to include 
the effect of non-Gaussianity of the cage shape in the theory. A first step has been performed in [150], 
but the detailed predictions of the non-Gaussian theory around jamming have not been worked out yet. 
Besides this, there are reasons to suspect that the exact description of the jamming transition within 
the replica theory requires an infinite number of replica symmetry breakings [76]. First of all, the 
generic situation with systems which have a glass transition to a one-step replica symmetry breaking 
(1RSB) solution describing the glass is to have a further transition to a phase with more - eventually 
infinite - levels of replica symmetry breaking. In mean field p-spin glasses, IRSB equilibrium states 
have a spectrum with no soft modes - and this is true of all but the very highest metastable states. 
A full replica symmetry breaking scheme would instead naturally bring in soft modes, as happens for 
example in the case of the Sherrington-Kirkpatrick model (the analogy between the physics of soft 
modes and those of the SK model has been already noted in [171]). Performing an infinite replica 
symmetry breaking calculation in the context of particle systems seems a difficult task: however, the 
calculation might be doable in the limit of d — > oo along the lines of [150]. A detailed analysis of the 
Bethe lattice model of [66] should help to shed light on this difficult problem. Also, setting up a field 
theory description of the jamming transition, using the replica method, along the lines of [124], might 
be possible. 
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